AD-A250  801 


NUMERICAL  STUDY  OF  AXIAL  TURBULENT  FLOW 
OVER  LONG  CYLINDERS 


Joao  C.  Neves.  Parviz  Moin  and  Robert  D.  Moser 


Prepared  with  the  support  of  the 
Office  of  Naval  Research 
under  contract  N00014-88-k-0145 
and  in  cooperation  with  the 
Naval  Underwater  Systems  Center 


s 


DTIC 

ELECTE  I 
JUU011392  I 


Report  No.  TF-54 


Thermosciences  Division 
Department  of  Mechanical  Engineering 
Stanford  University 
Stanford.  California  94305 


April  1992 


!  92  137 


NUMERICAL  STUDY  OF  AXIAL  TURBULENT  FLOW 


OVER  LONG  CYLINDERS 


By 


Joao  C.  Neves,  Parviz  Moin  and  Robert  D.  Moser 


Prepared  with  the  support  of  the 
Office  of  Naval  Research 
under  contract  N00014-88-k-0145 
and  in  cooperation  with  the 
Naval  Underwater  Systems  Center 


Report  No.  TF-54 


Thermosciences  Division 
Department  of  Mechanical  Engineering 
Stanford  University 
Stanford,  California  94305 


Accesion  For 


NTIS  CRA&I 
DTiC  TAB 
UnaiKiouDced 


Justdication 


By 


Distribution/ 


AvaHability  Codes 

Dist 

fi-l 

Avail  < 
Spe 

i!.d,'or 

cial 

April  1992 


Statement  A  per  telecon  James  Fein 

ONR/Code  1221 

Arlington,  VA  22217-5000 


NWW  5/29/92 


/ 


AD-A250801 


Woids/Phrases(4 

TEXT 


v/ords  max)  that  match  Tii-^sauiu: 
THESAURUS 


Ell  tb.  ie  s 


@1@ 

§200 

0230 


@240 

@270 


AD-A250801 

u 

Boundaries 
Equations , 
Layers,  He< 
Number,  Sci 

Walls , 
u 

Convex  trai 
significan 


BOUNDARY 
BOUNDARY  LAYER 
CABLES 
CHANNEL 
CURVATURE 
ENERGY 
EQUATION 
FLOW 
FLOW  FIELDS 
FRICTION 
KINETIC 
KINETIC  ENERGY 
LAYER 
LAYERS 
MEAN 
MOMENTUM 
NUMBER 
PRESSURE 
QUADRANT 
RE( 

REGIONS 
REYNOLDS  NUMBER 
SCALE 
SIMULATION 
SIMULATIONS 
SKIN  FRICTION 
SLOPE 
SONAR 
SPEED 
STATISTICS 
THICKNESS 
TRANSVERSE 
TURBULENCE 
TURBULENT  FLOW 
TURBULENT  FLOWS 
VELOCITY 
VORTICES 
WALL 

92-17 


Boundaries 

Boundary  Layer 

Cables 

Channels 

Curvature 

Energy 

Equations 

Flow 

F 1  o  v;  Fields 

Friction 

Kinetics 

Kinetic  Energy 

Layers 

Layers 

Mean 

Momentum 

Numbers 

Pressure 

Quadrants 

Regions 

Regions 

Reynolds  Number 
Scale 

Simulation 
Simulation 
Skin  Friction 
Slope 
Sonar 
Velocity 
Statistics 
Thickness 
Transverse 
Turbulence 
Turbulent  Flow 
JOTurbulent  Flow  j 
Velocity 
Vortices 
Walls 


\  /T'-'*  .7  Lty  5  S 


,  Boundary  Layer,  Channels,  Curvature,  Energy, 

Flow  Fields,  Fi  ic  ti  ine  t  i  c  Energy  ,  Kin-et.ic»-er^ 

an.  Momentum,  lUimbe Pressure,  Quadrants,  Regions,  Reynolds 
ale,  S-imul at Skin  Friction,  Slope,  Sonar,  Statistics, 
Transverse,  Turbulence  ,}^urbulent  Flow,  Velocity,  Vortices, 

nsverse  curvature  effects  in  wall  bounded  turbulent  flows  are 
t  if  the  boundary  layer  thickness  is  large  compared  to  the 


AD-A250801  (Cont'd) 


radius  of  curvature.  The  curvature  affects  the  inner  part  ot  the  flow  if 
a+,  the  cylinder  radius  in  v/all  units,  is  small.  This  flow  regime  is 
common  in  sonar  devices  towed  by  long  cables.  Two  direct  numerical 
simulations  of  transversely  curved  flows  v;ere  performed  and  statistical 
and  structural  data  were  extracted  from  th=  computed  flow  fields.  The 
effects  of  the  transverse  curvature  were  identified  by  ccmparing  the 
present  results  v;ith  those  of  the  plane  channel  simulation,  performed  at 
a  similar  Reynolds  number.  As  expected,  the  transversely  curved 
turbulent  flow  exhibits  many  features  common  to  the  planar  flov;&: 
near-wall  lov;  speed  streaks.  near-v;all  inclined  shear  layers,  near-v;all 
streamwise  vortices,  Reynolds  shear  stress  dominated  by  second  and 
fourth  quadrant  events,  etc.  As  the  curvature  increases,  the  skin 
friction  increases,  the  slope  of  the  logarithmic  region  decreases  and 
turbulence  intensities  are  reduced.  Several  turbulence  statistics  are 
found  to  scs.]'*  *’ith  a  curvature  dependent  velocity  scale  derived  from 
the  mean  momentum  equation.  Near  the  vjall,  the  flov;  is  more  anisotropic 
than  in  the  plane  channel  with  a  larger  percentage  of  the  turbulent 
kinetic  energy  resulting  from  the  streamwise  velocity  fluctuations. 
Near-wall  streamwise  vortices  are  the  strongest  sources  of  pressure 
fluctuations.  As  the  curvature  increases,  regions  of  strong  normal 
vorticity  develop  near  the  wall. 
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Abstract 


Convex  transverse  curvature  effects  in  wall  bounded  turbulent  flows  are  sig¬ 
nificant  if  the  boundary  layer  thickness  is  large  compared  to  the  radius  of  curvature. 
The  curvature  affects  the  inner  part  of  the  flow  if  a"*",  the  cylinder  radius  in  wall 
units,  is  small.  This  flow  regime  is  common  in  sonar  devices  towed  by  long  cables. 

Two  direct  numerical  simulations  of  transversely  curved  flows  were  performed 
and  statistical  and  structural  data  were  extracted  from  the  computed  flow  fields. 
The  effects  of  the  transverse  curvature  were  identified  by  comparing  the  present 
results  with  those  of  the  plane  channel  simulation,  performed  at  a  similar  Reynolds 
number.  As  expected,  the  transversely  curved  turbulent  flow  exhibits  many  fea¬ 
tures  common  to  the  planar  flows:  near-wall  low  speed  streaks,  near-wall  inclined 
shear  layers,  near- wall  streamwise  vortices,  Reynolds  shear  stress  dominated  by  sec¬ 
ond  and  fourth  quadrant  events,  etc.  As  the  curvature  increases,  the  skin  friction 
increases,  the  slope  of  the  logarithmic  region  decreases  and  turbulence  intensities 
are  reduced.  Several  turbulence  statistics  are  found  to  scale  with  a  curvature  de¬ 
pendent  velocity  scale  derived  from  the  mean  momentum  equation.  Near  the  wall, 
the  flow  is  more  anisotropic  than  in  the  plane  channel  with  a  larger  percentage 
of  the  turbulent  kinetic  energy  resulting  from  the  streamwise  velocity  fluctuations. 
Near-wall  streamwise  vortices  are  the  strongest  sources  of  pressure  fluctuations.  As 
the  curvature  increases,  regions  of  strong  normal  vorticity  develop  near  the  wall. 

As  the  curvature  increases  the  wall  pressure  fluctuations  are  dominated  by 
smaller  spanwise  length  scales.  Fractional  contributions  from  various  layers  in  the 
flow  to  the  wall  r.m.s.  pressure  are  marginally  affected  by  the  curvature.  The  ratio 
of  the  axial  to  azimuthal  length  scales  of  the  wall  pressure  fluctuations  increase  with 
increasing  curvature.  Curvature  dependent  time  and  length  scales  are  identified 
that  collapse  the  high  frequency  range  of  the  wall  pressure  temporal  spectra  and 
the  high  wave  number  range  of  the  wall  pressure  streamwise  spectra  of  flows  with 
different  curvatures.  Taylor’s  hypothesis  holds  for  the  wall  pressure  fluctuations 
with  a  lower  convection  velocity  than  in  the  planar  case. 
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CHAPTER  1 


Introduction 


Turbulent  flows  that  evolve  over  surfaces  with  convex  curvature  normal  to  the 
mean  flow  are  common  in  engineering  applications.  However,  it  is  only  for  strong 
curvatures  that  the  curvature  effects  become  noticeable.  For  this  reason  this  flow 
has  received  less  attention  than  its  planar  counterpart,  and  the  body  of  experimen¬ 
tal  data  available  is  limited.  Transversely  curved  turbulent  flows  with  large  curva¬ 
ture  effects  occur  over  sonar  devices  towed  by  long  cables.  Of  particular  interest 
is  the  characteristics  of  the  wall  pressure  fluctuations  underneath  the  transversely 
curved  boundary  layers. 

1.1  Transverse  Curvature  Effects 

In  turbulent  flows  with  transverse  curvature  there  is  an  additional  length  scale, 
the  cylinder  radius,  a.  The  added  complexity  introduced  by  the  new  length  scale 
is  evident  in  the  laminar  flow  regime.  While  in  the  planar  case  the  absence  of  a 
length  scale  leads  to  the  self-similar  Blasius  velocity  profile,  in  the  axisymmetric 
boundary  layer  there  is  no  self-similar  laminar  solution  (Seban  &  Bond  [1951]). 

The  new  length  scale  gives  rise  to  several  flow  regimes  that  are  characterized 
by  the  ratios  of  the  cylinder  radius  to  the  flow  length  scales:  the  boundary  layer 
thickness  6,  and  the  viscous  length  scale  iz/ut-  The  two  resulting  parameters,  S/a 
and  a'*',  define  a  two-dimensional  parameter  space  in  which  three  flow  regimes  can 
be  identified.  If  S/a  is  small  the  curved  boundary  layer  is  similar  tc  the  planar 
boundary  layer.  In  some  applications  the  boundary  layer  may  grow  to  be  several 
times  thicker  than  the  radius  (large  S/a),  in  which  case  the  transverse  curvature 
affects  the  flow  differently  depending  on  the  magnitude  of  a'*'.  If  6/a  is  large,  and 
a'^  is  also  large,  then  the  curvature  only  affects  the  outer  part  of  the  flow.  If  S/a 
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Figure  1.1  Experimental  data  available  for  the  axial  turbulent  flow  along  a  circular 
cylinder  given  as  a  function  of  the  parameters  S/a  and  :  a  Luxton  et  al.  (1984); 
+  Lueptow  et  al.  (1985)  and  Lueptow  &  Haritonidis  (1987);  x  Willmarth  «  Yang 
(1970)  and  Willmarth  et  al.  (1976);  o  Rao  &  Keshavan  (1972);  o  Afzal  &  Narasimha 
(1976);  V  Richmond  (1957);  n  Yu  (1959);  •  present  calculations. 

is  large  and  a"*"  is  small,  the  curvature  affects  both  the  inner  and  the  outer  parts  of 
the  flow. 

Some  of  the  flow  parameters  6 /a  and  a'*'  that  have  been  investigated  experi¬ 
mentally  axe  summarized  in  Figure  1.1.  Lueptow  (1988)  provides  a  comprehensive 
review  of  the  experimental  investigations  of  this  flow.  With  the  exception  of  Afzal  & 
Narasimha  (1976),  whose  boundary  layers  have  characteristics  very  similar  to  those 
of  the  flat  plate,  in  most  experiments  the  influence  of  curvature  is  predominantly 
on  the  outer  part  of  the  flow.  Some  of  the  experiments  of  Luxton  et  al.  (1984) 
(9  <  a'*'  <  47  and  26  <  S/a  <  42)  2uid  of  Willmarth  et  al.  (1976)  {2  <  S/a  <42  and 
>  33)  are  within  the  range  where  both  the  inner  8ind  outer  parts  of  the  flow  are 
affected.  In  the  tAperiments  of  Rao  &:  Keshavan  (1972)  (4  <  6/a  <  12  and  o'*"  >  23) 
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and  Lueptow  et  al.  (1985)  and  Lueptow  &  Haritonidis  (1987)  (4  <  6/a  <  8  and 
>  38)  the  curvature  primarily  affects  the  outer  layer. 

One  of  the  characteristics  of  flows  with  transverse  curvature  is  larger  skin  friction 
coefficients  (Cf)  than  in  planar  flows  of  similar  Reynolds  number  (see  for  example 
Rao  &  Keshavan  [1972]  or  Willmarth  et  al.  [1976]).  In  the  viscous  sublayer  of 
transversely  curved  flows  the  momentum  equilibrium  is  expressed  by  rr  =  ctw 
(Glauert  &  Lighthill  [1955]).  Reid  &  Wilson  (1963)  and  Rao  (1967)  proposed  a 
curvature  dependent  law  of  the  wall, 

=a+ln^l  +  ^^  ,  (1.1.1) 

which  predicts  curvature  effects  on  the  sublayer  for  sufficiently  small  . 

To  achieve  large  6 /a,  experiments  are  typically  performed  over  long  tubes  or 
wires  of  very  small  diameter  (Luxton  et  al.  [1984]  and  Lueptow  et  al.  [1985] 
and  Lueptow  Sz  Haritonidis  [1987]).  In  such  experimental  facilities,  the  structural 
isolation  of  the  wires,  the  aeroelastic  interaction  between  the  flow  and  the  wire, 
the  alignment  of  the  cylinder  with  the  mean  flow,  and  the  cylinder  sag  are  major 
concerns.  In  addition,  if  a'*'  is  small,  the  size  of  the  measuring  probe  {e.g.,  hot 
wire)  relative  to  the  cylinder  diameter  becomes  an  issue  in  near- wall  measurements. 
Finally,  some  of  the  measurement  techniques  used  are  based  on  the  assumption  that 
close  to  the  wall  the  mean  velocity  profile  is  the  same  as  that  of  the  flat  plate  (see 
for  example  Richmond  [1957]  and  Lueptow  et  al.  [1985]).  Attempts  to  address 
some  of  these  problems  are  described  by  Willmarth  &  Yang  (1970)  and  Lueptow 
&  Haritonidis  (1987). 

If  6/a  is  not  large  (Rao  &  Keshavan  [1972]  and  Lueptow  et  al.  [1985]),  the 
logarithmic  region  of  the  mean  velocity  profile  has  the  same  slope  as  the  planar 
case.  This  flow  regime  is  described  analytically  by  Afzal  &  Narasimha  (1976) 
with  the  method  of  matched  asymptotic  expansions.  They  conclude  that  when 
6 /a  =  0(1)  and  in  the  limit  of  a'*’  oo,  the  axisymmetric  turbulent  boundary 
layer  must  have  logarithmic  and  velocity  defect  profiles  of  the  same  form  as  in  the 
planar  turbulent  boundary  layer.  The  two  parameters  in  the  logarithmic  velocity 
profile  are  found  to  depend  on  6/a  (Afzal  &  Narasimha  [1976]).  As  6/a  increases 
and  o'*"  decreases,  the  mean  velocity  profiles  exhibit  a  logarithmic  region  with  a 
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decreasing  slope  (Lueptow  et  al.  [1985]),  and  for  sufficiently  large  6/a  and  small 

,  the  logarithmic  region  deteriorates  and  becomes  negatively  curved  ( Willmaxth 
et  al.  [1976]  and  Luxton  et  al.  [1984]).  Some  effort  has  been  devoted  to  identifying 
the  curvature  dependence  of  the  logarithmic  velocity  profile,  however,  it  is  clear 
from  Leuptow’s  (1988)  review  that  no  consensus  exists  on  this  issue.  Nevertheless, 
the  available  data  suggests  that  when  6/a  is  small  and  a'*'  is  large  the  slope  of  the 
logarithmic  region  tends  to  scale  with  6/a  (Lueptow  et  al.  [1985]).  As  6/a  increases 
and  a^  decreases,  the  logarithmic  profile  depends  on  both  parameters  (Willmarth 
et  al.  [1976],  Luxton  et  al.  [1984]). 

The  Reynolds  shear  stress  in  the  outer  part  of  the  transversely  curved  boundary 
layer  is  lower  than  its  planar  counterpart  (Lueptow  et  al.  [1985]).  Reynolds  shear 
stress  quadrant  contributions  in  the  transversely  curved  turbulent  boundary  layers 
are  similar  to  those  in  the  planar  geometry  (Lueptow  &  Haritonidis  [1987]).  At 
y'^  «  39  Lueptow  et  al.  [1987]  measured  higher  fractional  contributions  from  second 
and  forth  quadrant  events  than  in  the  planar  case  for  weak  events.  Intense  events 
arf'  found  to  have  a  lower  contribution  to  the  Reynolds  shear  stress. 

The  axial  and  normal  turbulence  intensities  measured  by  Luxton  et  al.  (1984) 
or  Lueptow  &:  Haritonidis  (1987)  are  also  lower  than  their  flat  plate  counterparts 
in  the  outer  part  of  the  boundary  layer.  Close  to  the  wall  {y'*'  <  20),  the  measured 
axial  intensities  have  magnitudes  similar  to  those  of  the  planar  case,  with  maxima 
/ut  ^  3.2)  also  located  at  i/'^  ~  12. 

In  a  visualization  study  of  axial  flow  over  a  cylinder,  Lueptow  &  Haritonidis 
(1987)  observed  large  scale  structures  moving  across  the  cylinder.  It  is  not  clear 
whether  the  observed  transverse  flow  motions  are  the  result  of  oscillations  of  the 
wire  (Lueptow  [1989]).  Since  these  large  structures  were  linked  to  the  azimuthal 
velocity  component  (Lueptow  &  Haritonidis  [1987])  they  should  significeintly  ziffect 
the  intensity  of  the  azimuthal  velocity  fluctuations,  which  unfortunately  were  not 
measured.  Simil2u:  large  scale  structures  were  zdso  observed  by  Luxton  et  al.  (1984) 
for  larger  curvatures  (6 /a  >  20)  only.  Luxton  et  al.  (1984)  suggested  that  these 
large  scale  structures  are  important  for  turbulence  generation  in  this  flow. 
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1.2  The  Wall  Pressure 


Knowledge  of  the  characteristics  of  the  wall  pressure  fluctuations  is  important 
for  imderstanding  flow-induced  sound  generation  and  the  structural  interaction 
between  the  flow  zind  the  body.  A  comprehensive  treatment  of  this  subject  is 
given  by  Blahe  (1986).  Since  transversely  cmrved  turbulent  flows  occur  over  sonar 
devices  towed  by  long  cables,  it  is  important  to  understand  the  characteristics 
of  the  wall  pressure  fluctuations  in  the  flow  regime  where  8 /a  is  large  and  a'*'  is 
small.  Reviews  of  measured  wall  pressure  fluctuations  on  a  flat  plate  are  given  by 
Willmarth  (1975)  and  Eckelmann  (1989).  The  effects  of  transverse  curvature  on  the 
space-time  characteristics  of  the  wall  pressure  were  investigated  experimentally  by 
Willmarth  &  Yang  (1970)  for  6/a  w  2  and  by  Willmarth  et  al.  (1976)  for  6/a  w  4. 
Since  in  both  studies  a"*"  =  0(10^),  the  curvature  effects  were  limited  to  the  outer 
part  of  the  flow. 

Many  of  the  difficulties  in  the  measurement  of  the  wall  pressure  fluctuations  can 
be  overcome  by  the  use  of  direct  numerical  simulations  (Handler  et  al.  [1984]  and 
Choi  &  Moin  [1990]).  However,  simulations  are  limited  to  low  Reynolds  numbers. 
There  is  some  evidence  (Willmarth  [1975])  that  the  root-mean-square  (r.m.s.)  pres¬ 
sure  normalized  by  the  mean  wall  shear  decreases  with  decreasing  Reynolds  number. 
This  was  confirmed  by  Choi  &:  Moin  (1990),  who  also  found  that  the  spectrum  of 
the  wall  pressure  fluctuations  decreases  with  decreasing  Reynolds  number  for  low 
frequencies  when  scaled  with  inner  variables  and  for  high  frequencies  when  scaled 
with  the  outer  variables.  In  the  two  transversely  curved  turbulent  boundary  layer 
studies  of  the  wall  pressure  fluctuations  no  appreciable  effect  of  curvature  on  the 
r.m.s.  wall  pressure  was  found  (Willmarth  &  Yang  [1970]). 

As  had  been  found  in  previous  experiments  (Willmarth  &  Wooldridge  [1962], 
Wills  [1964],  Emmerling  [1974]),  the  large  scale  pressure  correlation  contours  of 
Choi  &  Moin  (1990)  were  more  elongated  in  the  spanwise  direction  than  in  the 
streamwise  direction.  The  curvature  dependence  of  the  wall  pressure  correlations 
has  not  been  established.  Willmarth  &  Yang  (1970)  (6/a  w  2)  suggest  that  for  large 
separations  the  ratio  of  spanwise  to  streamwise  length  scales  decreases.  However, 
their  later  study  for  6/a  4  (Willmarth  et  al.  [1976])  shows  the  opposite  trend. 

Nevertheless,  the  wall  pressure  (Willmarth  et  al.  [1976])  is  significantly  correlated 
around  the  cylinder  and  the  azimuthal  extent  of  the  correlation  increases  with 
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increasing  curvature.  Willmarth  &  Yang  (1970)  argued  that  the  main  effect  of 
the  transverse  curvature  is  to  reduce  the  spanwise  length  scale  of  the  large  eddies, 
located  in  the  outer  part  of  the  boundary  layer. 

In  the  planar  turbulent  boundary  layer  (Willmarth  &  Wooldridge  [1962])  as  well 
as  in  the  transversely  curved  turbulent  boundary  layer  (Willmarth  &  Yang  [1970] 
and  Willmarth  et  al.  [1976]),  the  space-time  wall  pressure  correlations  indicate  that 
the  eddies  decay  after  traveling  a  distance  of  the  order  of  their  size.  Convection 
velocities  can  be  defined  from  the  wall  pressure  space- time  correlation  contours 
(Wills  [1964]).  In  the  planar  geometry,  there  is  a  consensus  among  the  experimen¬ 
tal  (Willmarth  &  Wooldridge  [1962],  Panton  &  Linebarger  [1974])  and  numerical 
(Choi  &  Moin  [1990])  values  reported  for  the  convection  velocity,  which  is  about 
O.Sf/oo-  Using  a  definition  of  the  convection  velocity  that  is  independent  of  tempo¬ 
ral  separation  (Wills  [1964]),  Choi  &  Moin  (1990)  found  a  convection  velocity  of 
0.72I7oo,  for  which  Taylor’s  hypothesis  was  more  accurate  for  large  structures  than 
for  small  ones.  In  the  transversely  curved  turbulent  flows  of  Willmarth  &  Yang 
(1970)  and  Willmarth  et  al.  (1976),  convection  velocities  identical  to  that  of  the 
planzu*  case  were  reported. 

1.3  Motivation  and  Objectives 

Direct  numerical  simulation  of  turbulence  is  the  ideal  tool  for  the  study  of  turbu¬ 
lent  flows  at  low  Reynolds  numbers.  This  technique  is  also  well  suited  for  providing 
spatially  and  temporedly  resolved  data  to  study  the  characteristics  of  the  wall  pres¬ 
sure  fluctuations,  as  well  as  for  identifying  the  structiural  features  in  the  flow. 

Since,  as  the  curvature  increeises,  the  circumference  of  the  cylinder  decreases, 
there  will  be  fewer  structures  around  the  cylinder.  In  this  sense  the  transversely 
curved  flow  in  the  large  6 /a  and  small  regime  is  a  better  setting  for  the  study 
of  the  mechanisms  of  waJl  bounded  flows.  This  is  similar  to  the  minimal  channel 
of  Jimenez  &:  Moin  (91)  with  the  added  advantage  that  the  periodic  boundary 
condition  in  the  transverse  direction  axe  the  natural  boundary  conditions. 

The  study  described  in  the  following  chapters  had  the  following  objectives: 

»)  to  develop  a  pseudo-spectral  method  for  the  incompressible  Navier- 
Stokes  equations  in  cylindrical  coordinates  with  a  no-stress  outer  boundary; 
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«)  to  generate  a  spatially  resolved  database  for  turbulent  flow  with 
transverse  curvature  to  study  curvature  effects  on  the  turbulence  statistics 
as  well  as  instzintaneous  flow  structures; 

t«)  to  identify  scaling  parameters  that  link  the  statistics  of  the  trans¬ 
versely  curved  flows  with  those  of  the  planar  case; 

iv)  to  study  the  effects  of  the  transverse  curvature  on  the  pressure 
field;  to  generate  a  temporally  and  spatially  resolved  database  of  the  wall 
pressure  and  wall  shear  stresses  and  to  study  the  space-time  characteristics 
of  the  wall  pressure  fluctuations. 
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CHAPTER  2 


The  Model  Problem 


This  chapter  describes  the  numerical  approach  used  to  solve  the  incompressible 
Navier-Stokes  equations  for  the  axial  flow  over  a  cylinder.  The  key  features  of 
the  method  are  the  truncation  of  the  semi-infinite  domain  in  the  radial  direction, 
the  decoupling  of  the  viscous  operators  of  the  radial  and  azimuth^d  momentum 
equations,  the  calculation  of  the  pressvire  through  a  variant  of  the  method  of  Kleiser 
&  Schumann  (1981)  and  the  imposition  of  the  far  field  boundary  conditions.  In 
addition,  a  new  dependent  variable  is  defined  (T),  and  the  velocity  is  obtained 
from  the  solution  of  two  partial  differential  equations  for  Vr  and  T  and  algebraic 
equations  for  Vg  and  Vz.  The  validation  of  the  numerical  approach  is  discussed  at 
the  end  of  the  chapter  with  the  calculation  of  the  stability  characteristics  of  several 
two-  and  three-dimensional  Taylor-Couette  flows. 

Throughout  the  chapter  the  velocity  scale  is  the  friction  velocity,  Ur  =  s/tw/p-, 
where  r«;  is  the  mean  wall  sheeir  stress.  The  thickness  of  the  layer  is  6,  and  the 
superscript  denotes  the  scaling  with  the  viscous  length,  vfur-  Total  fields  are 
indicated  by  upper  case  symbols  and  fluctuating  fields  are  denoted  by  lower  case 
symbols.  For  example,  V  and  v  are  the  total  and  fluctuating  velocity  vectors, 
respectively. 

2.1  Equations  of  Motion 

To  solve  the  Navier-Stokes  equations  for  axial  flow  along  a  cylinder  we  use  cylin¬ 
drical  coordinates  (r,  0, 0),  where  r  is  the  radial  coordinate,  6  is  the  azimuthal  co¬ 
ordinate  and  z  the  axial  coordinate.  In  cylindrical  coordinates  the  incompressible 
Navier-Stokes  equations  are  (see  for  example  Batchelor  [1967]); 


dVr 

dt 


+  Hr 


dP  J_ 
dr  ^  Rct 


( 


2  dVg\ 
d9  )' 


(2.1.1  a) 
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dt 


IdP  1 

r  do  ^  Rcr 

(2.1.1  h) 

dVz 

dz  + 

(2.1.1  c) 

where  H  =  {Hz,  Hr-,  He)  are  the  convective  terms  and  the  density  is  absorbed  into 
the  pressure.  Note  that,  unlike  their  Cartesian  counterparts,  the  viscous  operators 
of  the  radial  and  azimuthal  momentum  equations  are  coupled.  The  continuity 
equation  is  given  by 


^  1^ 

r  dr  r  do  dz 

The  Laplacian  operator  is 


(2.1.1  d) 


^  1  5^  d‘^ 

= - 1 - 1 - 1 - 

dr^  rdr  dO^  dz^' 

and  the  convective  terms  {Hz,  Hr, He)  are 


Hr  =  Vr 


dVr 

dr 


+ 


r  do  ^  ^  dz 


He 


..  dVe  ^  Ve  dVe 
^  dr  r  do 


dVe 


V.Vr 

r  ’ 


Hz 


^dVz  ^VedVz 
^  dr  r  do 


+  Vz 


dVz 

dz 


(2.1.2) 


(2.1.3) 


From  the  divergence  of  the  momentum  equations  one  obtains  the  Poisson  equa¬ 
tion  for  the  pressure 

v2p  =  _V.H.  (2.1.4) 

The  flow  under  consideration  is  naturally  homogeneous  in  the  azimuthal  direction 
and  is  assumed  to  be  homogeneous  in  the  axial  direction.  In  the  radial  direction  the 
semi-infinite  domain  is  truncated  to  a  finite  one,  where  r  E  [a,a  -|- 1].  At  r  =  a  -f  1, 
model  boundary  conditions  are  imposed.  The  flow  is  driven  by  a  mild  streamwise 
pressure  gradient  zmd  reaches  a  statistically  steady  state.  This  model  corresponds 


10 


to  a  boundary  layer  with  no  spatial  growth  and  is  therefore  an  approximation  to 
the  physical  case.  Apart  from  the  numerical  advantages,  the  lack  of  spatial  growth 
affords  a  better  statistical  sample  for  turbulence  correlations. 

Note  that  the  domain  truncation  of  this  model  prevents  the  occurrence  of  large 
scale  viscous/inviscid  interactions.  Experiments  (Luxton  et  al.  [1984],  Lueptow 
&  Haritonidis  [1987])  have  reported  that  such  structures  periodically  cross  the 
cylinder.  As  discussed  in  Chapter  1  the  origin  of  these  structures  may  be  due  to 
vibrations  of  the  cylinder  and  not  purely  a  fluid  dynamic  phenomenon.  In  any  case, 
our  model  is  incapable  of  capturing  such  structures. 


2.1.1  Boundary  Conditions 


On  the  surface  of  the  cylinder  (r  =  a)  no-slip  boundary  conditions  are  imposed. 


V 


=  0. 


(2.1.5  o) 


In  the  axial  direction,  a  finite  domain  of  length  Lz  is  chosen  and  periodic  boundary 
conditions  are  imposed.  In  the  azimuthal  direction  the  flow  is  natur^dly  periodic. 
In  a  semi-infinite  domain  the  far  field  (y  oo)  velocity  boundary  conditions  for 
the  axial  flow  along  the  cylinder  are  Vz  and  {Vr,Ve)  -+  0.  In  the  truncated 

domain  in  which  the  calculations  are  performed,  we  require  that  the  radi^ll  compo¬ 
nent  of  the  velocity  be  zero  at  r  =  a  -f  1.  The  boundary  conditions  for  the  axial 
and  azimuthal  velocity  components  are  obtained  by  requiring  that  the  surface  shear 
stress  be  zero  at  the  outer  edge  of  the  domain,  that  is,  <TrSlr=a-H  =  ‘^rz|r=a-|-l  =  0. 
Thus,  we  impose 


Vr 


r=a-H 


=  0, 


dVz 

dr 


=  0. 

r=a-t-l 


(2.1.5  b) 


The  computational  domain  is  illustrated  in  Figure  2.1. 
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2.1.2  Radial  Coordinate  Stretching 

In  the  radial  direction  Chebyshev  polynomial  expansions  (see  for  example  Got¬ 
tlieb  &  Orszag  [1967])  are  used.  These  expansions  promote  a  high  grid  resolution 
at  the  edges  of  the  domain  (r  =  a  and  r  =  a  -t- 1).  This  is  appropriate  for  the  inner 
boundary  where  no-slip  boundary  conditions  cire  imposed.  However,  at  the  outer 
edge  of  the  domain  fine  resolution  is  not  necessary  because  no  sizable  velocity  gra¬ 
dients  are  expected  there.  To  avoid  this  waste  of  resolution,  the  radial  coordinate 
mapping 

r  =  (2.1.5) 

is  used,  which  increases  the  resolution  close  to  the  cylinder  surface  at  the  expense  of 
the  resolution  at  the  outer  edge  of  the  computational  domain.  In  this  expression, 
the  Chebyshev  variable  ^  €  [—1,1]  ^  r  £  [a,  a  4-  1]  and  the  parameter  0  =  6 
was  used.  Because  the  axial  and  azimuthal  directions  are  periodic,  we  use  Fourier 
expansions  in  these  directions. 

2.2  Uncoupling  of  the  Viscous  Operators 

In  the  discussion  that  follows,  the  solution  procedure  is  discussed  without  refer¬ 
ence  to  the  radial  coordinate  mapping  (Equation  2.1.6)  because  its  inclusion  would 
only  increase  the  algebraic  complexity  of  the  equations.  Also  in  the  remainder  of 
the  chapter  the  Fourier  transformed  (in  6  and  z)  equations  will  be  considered. 

In  the  solution  of  the  incompressible  Navier-Stokes  equations  it  is  computation¬ 
ally  convenient  to  eliminate  the  pressure.  The  fourth  order  formulation  of  (Kim 
ei  al.  1987)  or  the  use  of  divergence-free  expansion  functions  for  the  velocity  field 
(Moser  et  al.  1983)  are  examples  of  how  this  can  be  achieved.  In  Cartesian  coor¬ 
dinates,  after  eliminating  the  pressure,  the  momentum  equations  are  only  coupled 
through  the  nonlinear  terms.  In  the  cylindrical  coordinate  system  the  coupling  is 
in  both  the  convective  terms  and  the  viscous  terms  of  the  radial  and  azimuthal 
momentum  equations.  In  order  to  decouple  the  radial  and  azimuthal  momentum 
equations  we  define  T  as  ("  denotes  the  Fourier  transform) 

t  =  ik.Vg  -  imVz,  (2.2.1) 
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and  ^ec^lst  the  momentum  equations  in  the  following  form: 


dVr  ,  -  dP  .  1 


where 


and 


^  =  k,m{-  -  r)P  +  -^V^t 

eft  r  Ker 


\  Vr  2  -  \ 


Ht  =  ikzHi  -  imffz  ~  (  -^  +  irn  —  Vr  |  , 

Rtr  \ 


-2  1  d  f  d\  nP  2 


(2.2.2  a) 


(2.2.2  b) 


(2.2.3) 


(2.2.4) 


is  the  Fourier  transformed  Laplacian  operator.  The  axial  wave-number  is 


kz  =  ^k, 


(2.2.5) 


and  the  azimuthal  wave- number  is  m. 

In  the  numerical  approach  adopted,  we  solve  first  for  Vr,  so  it  is  convenient  to 
express  the  boundary  conditions  for  T  ^ls  a  function  of  Vr-  At  the  surface  of  the 
cylinder  we  have  T(r  =  a)  =  0.  At  the  outer  edge  of  the  dommn  the  boundary 
conditions  on  the  velocity  (Equation  2.1.5  6)  and  continuity  give 


(m^-l-rA:^)^  -  klT 


=  —kzTTl 


r=:a-(-l 


dVr 

dr 


(2.2.6) 


r=a-|-l 
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The  computational  advantage  of  this  approach  is  that  in  Fourier  space  the  con¬ 
tinuity  equation  (Equation  2.1.1  d)  and  the  definition  of  T  (Equation  2.2.1)  deter¬ 
mine  an  algebraic  system  of  equations  for  Vz  and  as  a  function  of  Vr  and  T.  The 
solution  of  the  system  is 


Vz  = 


m2  -f-  rk\ 


m2  -f  rk\ 


m\Vr  + 


kz{Vr  +  r^]+mr 


(2.2.7) 


which  is  easier  to  solve  than  the  partial  differential  equations  that  they  replace. 


2.3  Time  Advancement 

An  implicit  second  order  Crank-Nicholson  scheme  is  used  to  time  advance  the 
viscous  terms  in  these  equations.  The  convective- like  terms,  Hr  and  Hri  are  ad¬ 
vanced  explicitly  with  a  second  order  Adams- Bashforth  scheme.  Note  that  the 
viscous  contributions  to  Hr  and  Hr  dio  not  involve  radial  derivatives,  which  could 
lead  to  numerical  stability  problems.  A  von  Neumann  stability  analysis  of  the  cross¬ 
terms  in  the  viscous  operators  shows  that  they  do  not  produce  a  severe  stability 
limitation  and  can  be  treated  explicitly. 

The  time  discretized  equations  for  Vr  (Equation  2.2.2  a)  and  T  (Equation  2.2.2  b) 
at  time  step  n  -h  1  are: 


At 


~  (3^”  -  +(1  +  ^v2)F,",  (2.3.1  a) 


=  0, 


r=a-|-l 


=  0, 
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(1  -  =  kzmAt{-  -  r)P^ 

zRct  r 

-f  (3h;- »?->)+ (i  +  ^v^)t», 
t"+M  =0, 


dr 


n+l 


(m^  +  rfci)-^  —  = —mk 

dr 

' r=a+l 


'  dr 


r=a+l 


(2.3.1  b) 


V^P"  =  -v  (2H  -jH  ), 


(2.3.1  c) 


where  the  boundary  conditions  for  the  pressure  (Equation  2.3.1  c)  enforce  the 
continmty  equation  at  the  boundaries  (see  Section  2.4). 

In  this  formulation,  two  elliptic  operators,  (1  —  A<V^/(2jRer))  and  V^,  need  to 
be  inverted  for  each  wave  number  pair  (m,k)  at  eeich  time  step. 

The  flow  variables  are  represented  by  the  following  expansions 


Nr  N,/2-l 

( E 

^  ^  1=0  m=-N 


Jm9 


m=-Nt/2 


JV,/2-l 

k=-N,/2 


(2.3.2) 


where  7/(^)  are  the  Chebyshev  polynomials  and  V/^jt(<)  and  Pimki^)  the  com¬ 
plex  spectral  coefficients  for  the  velocity  and  pressure,  respectively.  {Nr,  N0,Nz) 
denote  the  number  of  modes  in  the  radial,  azimuthal  and  axial  directions  respec¬ 
tively. 

With  the  radial  coordinate  stretching,  each  elliptic  operator  results  in  a  matrix 
with  12  diagonals  above  the  main  diagonal  and  8  below.  In  addition,  the  boundary 
conditions  fill  the  two  bottom  rows  of  each  matrix.  The  high  computational  cost 
of  the  inversion  of  matrices  with  this  bauid  structure  was  one  of  the  reasons  for 
the  adoption  of  a  constant  time  step.  Substantial  savings  in  computer  time  are 
achieved  by  the  precalculation  and  storage  of  the  LU  decompositions  of  the  matrices 
corresponding  to  the  operators  (1  —  AtV^/(2i2er))  and  V^. 
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The  const;  at  time  step  wzis  chosen  in  accordance  with  numerical  stability  con¬ 
straints.  The  CFL  number 


C  FL  =  TT  At 


/  ivy  ,  M 

\  Ar  rA9  Az ) 


(2.3.3) 


max 


was  required  to  be  less  than  0.6  throughout  the  calculations. 


2.4  Continuity  at  the  Boundaries 

The  role  of  the  pressure  is  to  ensure  a  divergence  free  velocity  field;  thus  the 
appropriate  boundary  conditions  for  the  pressure  at  r  =  a  and  r  =  a  -|-  1  must 
be  chosen  such  that  continuity  is  satisfied  everywhere  including  at  the  botmdaries. 
Kleiser  &  Schumann’s  (1981)  algorithm  for  the  enforcement  of  continuity  is  used 
here. 

We  start  by  writing  P  and  Vr  as 


pn+l  =pn+l  +a»+lp^^, 

V?+‘  =Vr"'^' +  a’i+'Vrk,  +  o^+'Vr^^, 


(2.4.1) 


where  and  0:2^*  are  constants  to  be  determined  at  each  time  step.  Pi  and 


n+1 


Vrj  are  the  particular  solutions  and  satisfy 


..  n  -  - 

V^F/*  =  -V  •  H  , 


F/* 


r=a,a-|-l 


=  0, 


(2.4.2) 


and 


At 


=  0, 


(2.4.3) 


=  0, 

r=a 


Vr?*' 


r=a-l-l 
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respectively. 

Pfiy  and  P/,2  are  obtained  from: 


Phi  =  Phi 

*  ^  #•— /t  * 


r=a 


r=a4-l 


=  0, 


(2.4.4) 


v2p.  =o>  =o>  Ph 


r=a 


r=a+l 


=  1. 


In  order  to  satisfy  Equations  2.3.1  a  and  2.4.3,  Fr/j,  and  Fr/,2  are  the  solutions  of 


Vr»,  =  0,  Vrk, 


r=a 


r=a+l 


=  0, 


(2.4.5) 


(1  -  vhv  j.  -  -a<2^ 

2Pcr^  ~  dr  ^ 


Vrhi  ~ 
•*  r=a 


r=a+l 


=  0. 


With  a  constant  time  step  At  (used  in  the  present  calculations),  Pf^^,  Pf^^,  Fr/j,  and 
Vrh2  need  to  be  calculated  once,  stored  and  used  throughout  the  calculations. 

The  constants  and  are  obtained  at  each  time  step  from  the  the  con¬ 

straint  imposed  by  the  continuity  equation  at  the  boundaries.  The  radial  derivative 
of  the  continuity  equation  and  the  velocity  boimdary  conditions,  given  by  Equations 
2.1.5  a  and  2.1.5  b,  leaid  to  the  following  compatibility  conditions  at  the  boundaries 
of  the  computationed  domain: 


dv;^+^ 

dr 


=  0, 

r=a 


(2.4.6) 
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The  expansion  of  Vr  according  to  Equation  2.4.1  and  the  above  compatibility 
conditions  lead  to  the  following  system: 


which  determines  and  With  these  constants,  the  pressure  and  radial 

velocity  fields  that  satisfy  the  continuity  equation  are  uniquely  determined. 


2.5  The  Mean  Flow 

The  mean  flow  (corresponding  to  m  =  0  euid  Ar  =  0)  in  the  axial  direction  is 
governed  by  the  axial  mean  momentum  equation 


where  the  over  bar  denotes  average  over  6  and  z  planes. 

The  flow  is  driven  by  a  mild  2ixial  pressure  gradient  so  that  it  can  reach  a 
statistical  steady  state.  In  the  cylindrical  geometry  the  mass  flux  in  the  axial 
direction  is  given  by 

a+1 

M  =  2n  J  Vzrdr.  (2.5.2) 

a 
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Integration  of  the  mean  axial  momentum  equation  in  the  radial  direction  gives 


dz 


1  +  2a 


a+1 


rdr  — 


2a 


dV, 


Rcril  +  2a)  dr 


r=a 


(2.5.3) 


which  determines  the  mean  axial  pressure  gradient  necessary  to  maintain  a  constant 
mass  flux. 

The  only  solution  for  Vr  that  satisfles  the  mean  continuity  equation  and  the 
boimdary  conditions  is 

Vr  =  0.  (2.5.4) 

Finally,  the  mean  flow  in  the  azimuthal  direction  (V#)  satisfies 


dV, 

dt 


Ve 


r=a 


=  0, 


(2.5.5) 


2.6  Code  Verification 

The  method  described  above  was  extensively  tested.  The  eigenvalues  of  the 
Stokes  operator  with  homogeneous  boimdary  conditions  were  computed.  As  ex¬ 
pected,  the  eigenvalues  are  real  iind  negative.  Furthermore,  they  were  verified 
to  be  in  close  agreement  with  the  eigenvalues  computed  by  Moser  (1988)  using  a 
different  numerical  method. 

In  addition  several  Taylor  Couette  flow  ce^es  were  computed.  The  extensive 
experimental  and  computational  data  available  on  Taylor  Couette  flow  makes  it  ein 
ideal  test  case  for  the  code.  First,  the  critical  Reynolds  number  (iicc)  for  the  onset 
of  Taylor  vortices  was  computed  and  compared  with  the  existing  data.  Comparisons 
with  analytical  (DiPrima  &  Eagles  [1977])  eind  computational  (Moser  et  al.  [1983]) 
results  for  two  gap  widths  are  shown  in  Table  2.1  o.  In  this  table  and  the  remainder 
of  the  chapter  Ri  and  Ro  are  the  radii  of  the  inner  and  outer  cylinder,  d  =  Rp  —  Ri  is 
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the  gap  width  and  A  is  the  axial  wave-number  of  the  vortices.  The  two-dimensional 
Taylor  Couette  flow  is  also  characterized  by  the  torque 


G  ~  ’htRtfi 


r=r, 


(2.6.1) 


that  must  be  applied  to  the  inner  cylinder  to  drive  the  flow.  Calculated  torque 
values  for  gap  widths  and  Rejmolds  numbers  for  which  there  are  both  experiments 
(Donnelly  &  Simon  [I960])  and  computations  (Moser  tt  al.  [1983])  are  reported  in 
Table  2.1  6.  In  all  cases  the  two-dimensional  caJculations  of  the  present  code  are 
in  excellent  agreement  with  both  the  experiments  and  previous  calculations. 


n/ro  =  0.5  Ti/ro  =  0.95 

A/d  =  3.976  A/d  =  4.018 

DiPrima  &  Eagles  (1977)  ^  184.99  68.19 

Moser  ei  al  (1983)  *  185.99  68.2 

present  184  68.2 


Table  2.1  a  Critical  Reynolds  number  {Rtc)  for  2-D  Taylor  vortices 


Donnelly  &  Simon  (1960)  ^ 
Moser  et  al  (1983) 
present 


fi/'f'o  =  0.5,  Rej  =  78.8 
A/d  =  3.976 

1.479  X  10^ 

1.487  X  10^ 

1.486  X  10^ 


n/'^'o  =  0.95,  Red  =  195 
A/d  =  4.018 

5.26  X  10^ 

5.42  X  10^ 

5.43  X  10^ 


Table  2.1  6  Torque  for  2-D  Taylor  vortices 


*  Computations;  f  Stability  analysis;  1  Experiments. 
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ri/r„  =  0.868 
Rej  =  458.1 


r,/ro  =  0.875 
Red  =  243.5 


King  et  al  (1984)  0.334411, •  t  ,  0.3347ftj  *  0.375711i  * 

present  0.334311, •  0.375911,- 

Table  2.1  c  Wave  speed  of  3-D  Taylor  vortices  (X/d  =  3.0,  m  =  6) 


To  test  the  full  three-dimensional  code,  the  wave  speeds  of  wavy  Taylor  vortices 
(Coles  [1965])  were  calculated  for  several  gap  Reynolds  numbers  (Red),  aocial  wave 
lengths  (A/d)  and  wave  numbers  (m).  The  results  of  the  present  code  are  compared 
in  Table  2.1  c  to  the  data  of  King  et  al.  (1984),  which  includes  experiments  and 
computations. 

2.6.1  Numerical  Resolution 

The  spatial  resolution  parameters  describing  the  two  turbulent  flow  simulations 
reported  here  are  shown  in  Table  2.2. 

Note  that  the  effective  azimuthal  resolution,  =  2irr/N0,  depends  on  the  radial 
location,  with  the  coarsest  resolution  at  the  outer  edge  of  the  domain  (r  =  a  -|-  1). 
Because  of  this,  the  simulations  are  very  well  resolved  in  the  azimuthal  direction, 
near  the  wall.  At  the  outer  edge,  the  resolution  appears  to  be  adequate  as  shown  by 
the  velocity  spectra  in  Figure  2.2.  In  Figure  2.3,  the  axial  spectra  of  the  velocity 
fluctuations  at  two  distances  from  the  wall  are  shown.  In  both  simulations  the 
axi2il  resolution  appears  to  be  adequate  with  no  energy  accumulation  at  the  high 
wave-numbers. 

The  length  of  the  computational  domain  is  set  to  be  sufficiently  large  to  ensure 
that  the  velocity  is  decorrelated  for  separations  larger  than  half  of  the  computa¬ 
tional  box,  so  that  the  imposed  periodic  boundary  conditions  will  have  minimal 
effects  on  the  turbulence.  Axial  two-point  correlations  (Figiire  2.4)  show  that  this 
is  roughly  the  case. 


*  Computations;  f  Stability  analysis;  f  Experiments. 
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6/a 

5 

11 

Lz 

Air 

6ir 

Lz^ 

2676 

4656 

L0.'*’  =  27ra'*’ 

268 

141 

=27r(a++6+) 

1608 

1692 

A+  =  Lz^/Nz 

14 

14 

A+  =  Le^^/N, 

4 

1 

K  =  itfo 

25 

13 

(Nz^Nr^N,) 

(192,96,65) 

(320,96,129) 

Table  2.2  Grid  resolution  parameters 


Figure  2.2  Azimuthal  (spanwise,  6Ak0  =  6/r)  spectra  of  the  velocity  fluctuations 
normalizedby  6  2md  Ur  for  6/a  =  5  (a)  andfor6/a  =  11(6): - Ev^Vz', - Ev^vf 
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Figure  2.3  Axial  spectra  of  the  velocity  fluctuations  normalized  by  6  and  Ur  for  ■ 

8/a  =  5  {a  and  6)  and  for  6/a  =  11  (c  and  d):  - Ev,v,\ - - E^gve-  | 

I 

"  I 

I 
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CHAPTER  3 


The  Velocity  and  Vorticity 


This  chapter  describes  the  major  statistical  characteristics  and  time  dependent 
structure  of  the  velocity  field.  All  the  statistics  presented  are  obtained  from  a 
sample  collected  after  the  calculations  had  reached  the  statistically  steady  state, 
and  are  calculated  using  spatial  averaging  over  the  domain  in  the  B  and  z  directions 
and  temporsd  averaging  that  covers  2.5  S/ut  time  units  in  the  SJa  =  b  case  and 

2.1  S/ut  time  units  in  the  Sfa  =  11  case.  Throughout  the  chapter  the  results  are 
compared  to  those  of  the  plane  channel  of  Kim  ti  al.  (1987). 

In  Section  3.1  the  properties  of  the  mean  fiow  are  studied  and  compared  with 
the  available  experimental  data.  In  Section  3.2  turbulence  intensities  and  the 
Re3niolds  shear  stress  are  shown  to  decrease  with  increasing  curvature.  In  addition, 
a  curvature  dependent  local  velocity  scale,  which  collapses  several  flow  statistics, 
is  derived  from  the  momentum  equations.  The  structure  of  the  Reynolds  shear 
stress  is  discussed  in  Section  3.3.  Vorticity  statistics  are  presented  in  Section  3.4. 
Finally,  the  instantaneous  structural  characteristics  of  the  flow  are  presented  via 
contour  plots  in  Section  3.5. 

3.1  Mean  Velocity  Statistics 

Mean  flow  parameters  for  the  present  transversely  curved  flow  simulations  as  well 
as  for  the  plane  channel  {S/a  =  0)  are  reported  in  Table  3.1.  Because  the  Reynolds 
numbers  of  the  three  simulations  are  similar,  the  differences  in  Table  3.1  are  due 
to  the  trsmsverse  curvature.  For  comparison.  Tables  3.2  o,  3.2  h  and  3.2  c  contain 
the  mean  flow  parameters  of  the  experiments  of  Willmarth  et  al.  (1976),  Luxton 
et  al.  (1984)  and  Leuptow  &  Haritonidis  (1987).  In  agreement  with  experimental 
observations,  the  skin  friction  coefficient. 
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Rtr 
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Rea 

— 
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6*/S 

0.141 

0.154 

0.152 

0*1 6 

0.087 

0.123 

0.131 

H 

1.62 

1.25 

1.15 

Table  3 . 1  Mean  flow  parameters 


increases  with  increasing  curvature  (by  as  much  as  63%  for  6 /a  =  11)  when  com¬ 
pared  with  the  plane  channel  flow  with  comparable  Reynolds  number.  However, 
in  the  simulations  the  magnitude  of  this  increase  is  larger  than  has  been  observed 
experimentally,  probably  because  of  the  much  higher  experimental  Reynolds  num¬ 
bers  (Rcj  =  0(10^)).  In  Table  3.1  the  boundary  layer  displacement  and  momentum 
thicknesses  reported  for  the  present  transversely  curved  flow  simulations  are  defined 
by  (Luxton  tt  al.  [1984]) 


a/6-^-l  _ 

{a/8  +  S*/6f  -{a/6f  =  2  J  ^1  -  rdr,  (3.1.2  a) 

a/s 


and 


a/6+l  _  _ 

(a/s  +  eysf  -  (a/sf  =  2  j  (0  (l  - 

a/s 


(3.1.2  6) 
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6/a 

4.7 

16 

37.5 

a+ 
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198 

46 

Cf 

3.04  X  10'^ 

4.18  X  10“^ 

7.84  X  10-3 

Res 

90380 

69280 

27600 

Rea 

19230 

4330 

736 

Table  3.2 

a  Experimental  mean 

flow  parameters  (Willmarth  et  al.  [1976]) 

6/a 

26 

26.9 

41.6 

a^ 

47.4 

32.1 

13 

Cf 

7.3  X  10-3 

0.01 

0.017 

Res 

20386 

12252 

5962 

Rea 

785 

455 

140 

6*/6 

0.184 

0.185 

0.187 

e*/6 

0.181 

0.182 

0.182 

H 

1.014 

1.02 

1.03 

Table  3.2 

h  Experimental  mean 

flow  parameters  (Luxton  et  al. 

[1984]) 

6/a 

6.74 

7.16 

8.0 

a^ 

288 

144 

72 

Cf 

3.5  X  10-3 

3.8  X  10-3 

4.1  X  10-3 

Res 

46330 

23700 

12800 

Rea 

6419 

3209 

1605 

Table  3.2  c  Experimental  mean  flow  parameters  (Lueptow  Sz  Haritonidis  [1987]). 
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Figure  3.1  Mean-velocity  profiles: - plane  channel  (Kim  et  al.  [198’^);  cylinders 

with - 6 fa  —  5  and - 6 fa  =  11; .  planar  law  of  the  wall,  U'^  =  y'*’,  and 

log  law,  C/+  =  2.51n(y‘*')  +  5.5. 

respectively.  Even  though  several  different  definitions  for  the  displacement  and 
momentum  thicknesses  in  transversely  curved  flows  have  been  proposed  in  the  lit¬ 
erature,  there  are  several  reasons  to  adopt  the  above  definitions.  The  expression  fcr 
the  displacement  thickness  (Equation  3.1.2  o)  is  consistent  with  the  usual  m£iss  flux 
displacement  zurgument,  and  the  momentum  thickness  defined  by  Equation  3.1.2  b 
guises  in  the  axisymmetric  integral  momentum  equation.  In  this  sense,  the  two 
definitions  have  the  same  physical  meaning  gis  their  planar  counterpgurts,  and  in 
the  flat  plate  limit  {S/a  0),  Elquations  3.1.2  reduce  to  their  planar  counterparts 
(multiplying  by  6/a  before  taking  the  limit  6/a  —*  0).  Luxton  et  al.  (1984)  mea¬ 
sured  6*/ 6  and  6*/ 6  according  to  the  above  expressions  and  find  higher  values  than 
in  the  simulations.  This  is  probably  a  result  of  the  larger  6/am  the  experiments. 
The  shape  factors  {H  =  6*/6*)  measured  by  Luxton  et  al.  (1984)  are  around  unity, 
whereas  the  shape  factor  in  the  simulations  is  larger  than  one  and  decrezises  with 
increasing  6 /a  due  to  an  increase  in  the  momentum  thickness. 
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3.1.1  The  Mean  Velocity  Profile 


The  mean  velocity  profile  of  the  two  transversely  curved  flows  (^/a  =  5  and  11) 
are  compared  to  that  in  the  plane  channel  flow  in  Figure  3.1.  Because  of  the  small 
values  of  a"*",  the  viscous  sublayer  mean  velocity  profiles  axe  slightly  affected  by 
the  curvature  of  the  wall.  In  contrast  most  experiments  (Rao  &  Keshavan  [1972], 
Willmarth  et  al.  [1976],  Lueptow  et  al.  [1985]),  find  no  effect  of  curvature  in  the 
mean  velocity  profile  for  y'*'  <  20.  Clearly,  in  many  of  these  experiments  a'*’  is  too 
large  for  there  to  be  a  perceptible  iimer  layer  effect  (for  example  w  140  and 
6/a  ful  in  Lueptow  &  Haritonidis  [1987]).  However,  even  in  experiments  in  which 
the  values  of  a'*'  are  close  to  or  below  those  of  the  simulations  (for  exeimple  Luxton 
ei  al.  [1984]  and  Willmarth  et  al.  [1976]),  the  measured  mean  velocity  profiles 
agree  with  the  planar  velocity  profile  for  y"*"  <  20.  This  disagreement  between 
the  experiments  and  the  simulations  may  be  due  to  measurement  difficulties  close 
to  the  wall  of  the  cylinder.  In  addition,  it  should  be  noted  that,  in  many  of  the 
mean  velocity  measurements,  the  estimation  of  the  friction  velocity  is  estimated 
by  assmning  that  the  mean  velocity  profile  of  the  transversely  curved  flow  is  not 
affected  by  the  curvature  near  the  wall.  The  mean  velocity  profiles  of  the  present 
numerical  experiments  suggest  that  this  assumption  is  not  correct,  at  least  in  the 
flow  regime  of  interest  here,  that  of  large  6/ a  and  small  a'*'. 

In  the  planar  boundary  layer,  the  assumption  that  the  total  stress  is  dominated 
by  the  viscous  stress  in  the  viscous  sublayer  leads  to  the  planar  law  of  the  wall,  U'^  = 
y'*'.  For  the  axisymmetric  boundary  layer  in  the  absence  of  a  pressure  gradient,  the 
mean  momentum  equilibrium  in  the  viscous  sublayer  is  expressed  by 


as  first  noted  by  Glauert  &:  Lighthill  (1955).  Assuming  that  the  Reynolds  shear 
stress  is  negligible  in  comparison  to  the  viscous  stress  leads  to  the  law  of  the  wall 
for  axisymmetric  flows  (Reid  &  Wilson  [1963],  Rao  [1967]), 


(3.1.4) 
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Figure  3.2  Total  stress  normalized  by  ur'.  analytical  for  cylinders  with - ^/a  =  5 

and - 5/a  =  ll;  symbols  for  computed  residts  with  o  for  ^/o  =  5  and  with  a  for 

S/a  =  11. 

Note  that  to  second  order  in  t/"*",  Equation  3.1.4  is  given  by 

Thus,  the  law  of  the  wall  for  the  viscous  subalyer  in  transversely  curved  flows  is 
equal  to  the  planar  viscous  sublayer  law  of  the  wall,  U'^  =  plus  a  curvature 
dependent  correction.  Clearly,  if  a'^  is  large  the  mean  velocity  proflle  near  the  wall 
does  not  deviate  appreciably  from  the  planar  case.  The  mean  velocity  profiles  given 
by  Equation  3.1.4  are  better  approximations  to  the  computed  near- wall  velocity 
profiles  than  the  planeir  law  of  the  wall. 

As  observed  in  several  experiments  (Lueptow  et  al.  [1985],  Luxton  et  al.  [1984] 
and  Willmarth  tt  al.  [1976]),  the  slope  of  the  mean  velocity  profile  in  the  logarith¬ 
mic  region  decreeises  with  increzising  curvature  and  the  profiles  become  negatively 


3'^ 


curved.  However,  the  dependence  of  the  slope  of  the  logarithmic  region  on  the  cur¬ 
vature  is  larger  in  the  simulations  than  in  experiments  (Willmarth  et  al.  [1976]). 
The  difference  seems  to  be  associated  with  the  flow  regimes  that  were  investigated 
by  the  various  authors.  Lueptow  ei  al.  (1985),  for  example,  show  that  the  slope 
of  the  logarithmic  region  does  not  change  for  small  6 fa.  This  is  not  surprising 
because  for  small  Sja  the  curvature  effects  are  limited  to  the  outer  part  of  the 
boundary  layer  in  high  Reynolds  number  flows.  Based  on  the  same  measurements, 
Lueptow  et  al.  (1985)  also  argued  that  the  slope  of  the  logarithmic  region  is  not 
a  function  of  a^.  The  latter  conclusion,  however,  is  questionable  because  in  those 
experiments  is  consistently  large.  These  observations  are  in  general  agreement 
with  the  asymptotic  analysis  of  the  axisymmetric  momentum  equations  in  the  Umit 
of  fi/a  =  0(1)  and  a'*'  — »  oo  (Afzal  &  Narahsima  [1976]).  To  lowest  order  in  this 
limit,  the  axisymmetric  boimdary  layer  has  both  a  logarithmic  region  eind  a  outer 
velocity  defect  law.  For  smaller  a"*"  the  parameters  of  the  velocity  profile  in  the  two 
regions  depends  on  both  6 /a  and  a'*'.  The  measurements  of  Luxton  et  al.  (1984) 
for  flows  with  larger  S/a  than  in  the  present  simulations  but  with  similar  show 
logarithmic  region  slopes  comparable  to  those  reported  here.  The  lack  of  consensus 
on  this  issue  is  stimmarized  in  the  review  of  Lueptow  (1988).  A  parametric  study 
is  necessary  before  a  conclusion  cjm  be  reached,  but  it  seems  that  the  slope  of  the 
logarithmic  region  is  a  function  of  both  and  6 /a  when  both  the  inner  and  outer 
layers  are  affected. 

3.2  Turbulence  Statistics 

The  statistical  steady  state  is  characterized  by  an  equilibrium  between  the  mean 
total  stress  and  the  applied  streamwise  pressure  gradient.  For  the  flows  considered 
here  this  is  expressed  by 


_  1  dV z  2  a  (o -b  6)^  — 

+  =  P«r— ; 


Re  dr 


r  {a  +  Sy  — 


(3,2.1) 


Figure  3.2,  shows  the  total  stress  (computed  and  analytical)  for  the  present  simula¬ 
tions.  Note  that  the  computed  zmd  analytical  curves  coincide  in  both  simulations, 
indicating  that  the  statistical  steady  state  has  been  achieved. 
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Figure  3.3  Root-mean-square  velocity  fluctuations  normalized  by  Ut'.  (a)  Axial  I 

intensity,  {h)  Normal  intensity,  (c)  Azimuthal  intensity; - plane  channel  (Kim  ■ 

tt  al.  [1987]);  cylinders  with - S/a  =  5  and - S/a  =  11. 
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Figure  3.4  The  energy  partition  parameter  K*  =  2v'i/{y^  +  v\)‘. - plane  channel 

(Lee  tt  al.  [1990]);  cylinders  with - 6/a  =  5  and - 6 /a  =  11. 

3.2.1  Turbulence  Intensities 

The  turbulence  intensities  normalized  with  the  friction  velocity  decrease  through¬ 
out  the  layer  as  the  curvature  increases  (Figure  3.3).  The  streamwise  component  is 
the  most  energetic  and  the  location  of  its  maximum  moves  slightly  towards  the  wall 
as  the  curvature  increases.  The  smedler  turbulent  kinetic  energy  in  the  transversely 
curved  flows  can  be  attributed  in  part  to  the  smaller  surface  area  over  which  vor- 
ticity  fluctuations  cam  be  generated  relative  to  the  volume  of  turbulent  supported. 
The  cylinder  surface  is  apparently  not  less  efiicient  as  a  source  of  turbulent  kinetic 
energy;  rather  it  has  to  supply  a  larger  volume.  This  geometric  difference  is  one  of 
the  reasons  for  the  reduced  intensity  levels  in  the  transversely  curved  flows.  As  the 
curvature  increases  the  turbulence  intensities  in  the  outer  part  of  the  flow  decrease 
in  agreement  with  severaJ  experiments  (Lueptow  et  al.  [1985],  Lueptow  &  Haritoni- 
dis  [1987]  and  Luxton  ei  al.  [1984]).  For  experiments  with  comparable  a'*',  Luxton 
ei  al.  [1984]  report  a  maodmum  vaJue  of  the  streamwise  intensity,  '  /Uoo  ^  0.16, 
similar  to  the  simulation  result  of  0.15.  However,  for  the  larger  experiments  of 
Lueptow  &  Hau-itonidis  (1987),  the  maximum  value  of  vl  /ur  of  3.3  is  higher  than 


Figure  3.5  Pressure  strain  term  of  the  streamwise  intensity  (uz^)  budget  equation 

normsJized  v  and  Ur:  cylinders  with - Sfa  =  5  and - 8ja  =  11; - plane 

channel  (Mansoiu-  et  al.  [1988]). 

in  the  flat  plate,  contrary  to  the  simulation  results.  This  may  be  due  to  the  fact 
that  in  their  flows  a'^  is  too  large  for  there  to  be  am  inner  layer  curvature  effect. 

Similau-ly,  the  normal  and  azimuthal  intensities  are  lower  than  in 

the  planar  case.  However,  the  reductions  in  the  normal  and  2izimuthal  intensities 
are  greater  thain  that  of  the  axial  intensity.  The  energy  partition  parameter  (Lee 
etal.  [1990]), 


vl  +  vj' 


(3.2.2) 


shown  in  Figure  3.4,  is  a  measure  of  the  relative  contributions  to  the  turbulent 
kinetic  energy  of  the  streamwise  turbulence  intensity  and  the  intensities  normal  to 
the  mean  flow.  The  attenuation  of  the  normal  eind  azimuthal  velocity  fluctuations 
is  strongest  for  y"*"  <  30  and  increases  with  curvature. 

Figure  3.4  suggests  that  the  trzmsfer  of  energy  from  the  streamwise  velocity  com¬ 
ponent  to  the  other  two  velocity  components  is  strongly  deunped  as  the  curvature 
increases.  In  the  budget  equations  of  the  intensities  (see  Appendix  A)  it  is  the 
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Figure  3.6  Reynolds  shear  stress  normalized  by  Ur’.  cylinders  with - 6/a  =  5 

and - 6/a  =  11; - plane  channel  (Kim  et  al.  [1987]). 

pressure  strain  terms  that  are  responsible  for  the  intercomponent  energy  transfer. 
As  shown  in  Figure  3.5,  as  the  curvature  increases,  the  pressure  strain  term  of  the 
budget  equation  for  the  streamwise  intensity  decreases  significantly  throughout  the 
layer. 


3.2.2  Reynolds  Shear  Stress 

The  Reynolds  shear  stress  (Figxure  3.6)  is  also  reduced  by  curvature.  The  loca¬ 
tion  of  the  maximum  of  the  Reynolds  shear  stress  profile  is  a  function  of  6/a  eind 
moves  towards  the  wall  with  incre€ising  curvature.  However,  this  does  not  alfect 
the  position  of  the  maximum  in  the  production  {—VzVrdVz/dy)  of  turbulent  kinetic 
energy  (!/■*■  ~  12,  see  Figture  3.10  b).  In  the  outer  part  of  the  layer  the  measure¬ 
ments  of  Lueptow  et  al.  (1985)  show  that  the  Reynolds  shear  stress  decreases  with 
increasing  curvature  in  agreement  with  the  simulation  results.  In  the  inner  layer, 
the  measurements  of  the  Reynolds  shear  stress  reported  by  Lueptow  et  al.  (1985) 
for  the  transversely  curved  flows  do  not  differ  appreciably  from  those  of  the  pla¬ 
nar  case.  In  contrast  with  the  present  simulations  the  near-wall  maximum  of  the 
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Figure  3.7  Correlation  coefficient:  cylinders  with - 6/a  =  5  and - 6 /a  =  11; 

- plane  channel  (Kim  et  al.  [1987]). 

Reynolds  shear  stress  (—vjv^/u^)  reported  by  Leuptow  et  al.  (1985)  is  about  0.8, 
similar  to  the  plane  chaimel  result  of  Kim  et  al.  (1987)  (Figure  3.6). 

The  velocity  correlation  coefficient,  shown  in  Figure  3.7,  suggests  that  there  axe 
important  differences  both  between  the  two  transversely  curved  flows  as  well  as 
between  them  and  the  plane  channel.  Close  to  the  wall  {y/6  <  0.2),  the  streamwise 
and  normal  velocity  fluctuations  are  increasingly  better  correlated  as  the  curvature 
increases.  Away  from  the  w?dl,  the  large  curvature  case  {6 /a  =  11)  shows  a  sig¬ 
nificant  reduction  in  the  correlation  coefficient.  This  may  be  an  indication  of  flow 
stabilization,  although  for  6/a  =  11  the  flow  is  a  self  sustaining  turbulent  flow  and 
reaches  statistically  steady  state. 
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Figure  3.8  Velocity  scaling  function  6/a): - plane  channel  (6/a  =  0); 

cylinders  with - 6fa  =  h  and - 6/a  =11. 

3.2.3  Velocity  Scaling 

As  a  function  of  the  distance  to  the  wall  the  total  stress  (Equation  3.2.1)  can  be 
rewritten  as 

In  this  form  it  can  be  easily  compared  with  its  planar  counterpart 

r  =  p«?(l-y/6),  (3.2.4) 

which  is  a  linear  function  of  y/6.  These  two  expressions  suggest  the  definition  of  a 
new  velocity  scale  for  the  transversely  curved  flows 

u  =  ur/.F(y/6;  6/ a),  (3.2.5) 

where 


39 


Note  that  in  the  planai  limit  (S/a  —*  0),  J-(y/S;  6/ a)  —*  1  and  the  planar  velocity 
scale,  u  =  «r,  is  recovered.  Also,  for  small  curvature  (S/a  <  1),  .F  1  throughout 
the  layer,  (Figure  3.8). 

When  scaled  with  u,  the  turbulence  intensities  and  Reynolds  shear  stress  profiles 
of  the  two  curved  flows  collapse  (Figure  3.9  and  Figure  3.10  a).  The  collapse  with 
the  plane  channel  stresses  is  also  satisfactory  except  for  the  transverse  turbulence 
intensity.  Turbulence  production  and  viscous  dissipation  apparently  do  not  scale 
with  u  near  the  wall  (Figure  3.10  h  and  3.10  c). 

3.3  Quadrant  Analysis 

The  correlation  coefficient  profiles  (Figure  3.6)  suggest  that  the  Reynolds  shear 
stress  producing  events  are  strongly  affected  by  curvature.  In  the  transversely 
curved  flows,  the  pzurtition  of  the  Reynolds  shear  stress  among  the  four  quadrants 
is  very  similar  to  that  of  the  plane  channel,  as  shown  in  Figure  3.11.  It  is  a 
characteristic  of  wall  bounded  flows  that  the  Reynolds  she2ir  stress  is  dominated 
by  the  second  quadrant  events  (vz  <  0  and  Vr  >  0)  in  the  outer  layer,  and  by  the 
fourth  quadrant  events  (vz  >  0  and  Vr  <  0)  in  the  inner  layer.  The  crossover  point 
between  the  domains  of  second  and  fourth  quadrant  events  occurs  at  the  same 
distamce  from  the  wjiU  (y"*"  ~  12)  as  in  the  pleine  channel. 

The  fractional  contributions  to  each  of  the  quadrants  of  the  Reynolds  shear  stress 
are  shown  in  Figures  3.12  and  3.13  for  S/a  =  5  and  S/a  =  11,  respectively.  At  a 
given  distance  from  the  wall,  the  fractional  contributions  of  each  of  the  quadrants 
to  the  Reynolds  shear  stress  in  the  transversely  curved  flows  are  similar  to  the 
corresponding  fractional  Reynolds  shear  stress  contributions  in  the  plane  channel 
(Kim  et  al.  [1987]).  This  invariance  of  the  fractional  contributions  with  curvature 
means  that  Reynolds  shear  stress  events  of  any  intensity  contribute  the  same  frac¬ 
tion  of  the  total  Reynolds  shear  stress  in  the  transversely  curved  flows  as  in  the 
plane  channel.  This  is  contraury  to  the  measurements  of  Lueptow  &  Haritonidis 
(1987),  who  found  that,  in  transversely  curved  flows,  low  intensity  Reynolds  shear 
stress  events  accounted  for  a  larger  percentage  of  the  total  Reynolds  shear  stress 
than  in  the  planair  case. 
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Figure  3.9  Root-mean-square  velocity  fluctuations  normalized  by  u:  (a)  Axial 

intensity,  (b)  Normal  intensity,  (c)  Azimuthal  intensity; - plane  channel  (Kim 

et  al.  [1987]);  cylinders  with - S/a  =  5  and - S/a  =  11. 


Figure  3.10  Reynolds  shear  stress  (a),  Production  (6)  euid  dissipation  (c)  of  tur¬ 
bulent  kinetic  energy  normalized  by  u  and  «:  cylinders  with -  8ja  =  b  and 

- 6/a  =  11; - plsune  channel  (Mansour  e<  al.  [1988]). 
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Figure  3.11  Reynolds  shear  stress  for  each  quadrant  normalized  by  the  mean 

Reynolds  shear  stress;  (a)  S/a  =  5,  (6)  S/a  =11: . first; - second; - third; 

- fourth  quadrant. 


Figure  3.12  Fractional  contribution  to  —VzVr  from  each  quadrant  as  a  function  of 

threshold  for  S/a  =  5; - first; - second; - third; .  fourth  quadrants: 

(a)  t/+  =  8  (6)  =  13.1  (c)  y+  =  51.1. 


I 

I 
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Figure  3.13  Fractional  contribution  to  —VzVr  from  each  quadrant  as  a  function  of 


threshold  for  6 
(a)  y+  =  8.9  (5' 


/a  =  11; 

i)  y+  =  1 


=  11: - first; 


13.2  (c)  y+  =  51.7. 


second; - third; .  fourth  quadrants: 


The  skewness  factors  of  the  velocity  fluctuations,  shown  in  Figure  3.14,  indicate 
that  there  is  a  strong  effect  of  the  curvature  on  the  streamwise  and  normal  veloc¬ 
ity  fluctuations.  As  expected  from  the  reflective  symmetry  of  the  Navier-Stokes 
equations,  the  skewness  of  Vg  is  nearly  zero  everywhere.  Deviations  from  zero  are 
attributed  to  an  inadequate  statistical  sample. 

Away  from  the  wall  {y'*'  >  40),  the  skewness  of  Vr  is  positive  and  increases 
slightly  with  increasing  curvature,  as  was  also  observed  by  Luxton  et  al.  (1984). 
In  this  region  there  axe  also  strong  negative  streamwise  velocity  fluctuations,  as 
denoted  by  the  negative  skewness  of  Vz,  which  decrecises  with  increasing  curvature. 

Increasing  the  curvature  does  not  significantly  affect  the  skewness  of  the  axial 
fluctuations  (vz)  close  to  the  wall  (t/+  <  20).  On  the  other  hand,  for  t/+  <  30, 
there  are  significant  differences  in  the  skewness  of  the  normal  velocity  fluctuations 
(vr)  between  the  three  flows.  While  in  the  plane  channel  for  5  <  j/"^  <  30  the 
skewness  of  Vr  is  negative,  in  the  transversely  curved  flows  the  region  of  negative 
skewness  diminishes  and,  for  S/a  =  11,  the  skewness  of  Vr  is  positive  throughout 
the  layer.  For  y'*'  <  12,  where  fourth  quadrant  events  dominate  the  Reynolds  stress, 
the  skewness  of  Vz  is  positive. 

The  flatness  profiles  of  the  velocity  fluctuations  are  shown  in  Figure  3.15.  Near 
the  wjJl  (y'*’  <  20),  the  flatness  of  Vr  and  Vg  decrease  sheirply  with  increasing 
curvature.  The  flatness  of  Vz  eilso  decreases  in  the  neax-wall  region  (j/"*"  <  5).  For 
5  <  y'*'  <  30  the  flatness  of  Vz  is  not  affected  by  the  curvature  and  away  from  the 
wall  (y^  >  30)  it  increases  with  increasing  curvature. 


3.4  The  Vorticity 


In  cylindrical  coordinates  the 
Batchelor  [1967]) 

flz 

fir 

fij 


vorticity  components  are  given  by  (see  for  example 

^1  /d(rVg)  _ 
r[  dr  de  )' 

^  (3.4.1) 

r  de  dz' 

dVr  dVz 
dz  dr 
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Figure  3.16  Root-mean- square  vorticity  fluctuations  normalized  by  i/  and  Ur:  (a) 

Axial  intensity,  (b)  Normal  intensity,  (c)  Azimuthal  intensity; - plane  channel 

(Kim  et  at.  [1987]);  cylinders  with - 6/a  =  5  and - S/a  =  11. 
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8/a 

0 

5 

11 

'  e 

16.7 

13.8 

13.5 

ytnax 

21.0 

19.7 

20.1 

2.31 

1.26 

0.83 

0.139 

0.092 

0.061 

iZep 

242.6 

110.0 

70.2 

Table  3.3  Streamwise  vortex  parameters 


The  vorticity  intensities  normalized  by  the  mean  wall  shear  stress,  shown  in 
Figure  3.16,  decrease  with  increasing  curvature.  Unlike  the  velocity  intensities, 
the  vorticity  intensities  do  not  collapse  in  the  outer  part  of  the  layer  when  scaled 
with  the  velocity  scale  u.  As  in  the  plane  channel,  the  axial  vorticity  intensity 
exhibits  a  near-wall  local  minimum  and  local  maximum.  Kim  ei  al.  (1987)  linked 
the  locations  and  intensity  of  these  extrema  to  the  average  position  and  strength  of 
the  near- wall  streamwise  vortices.  In  their  Rankine  vortex  model,  the  mean  radius 
is  estimated  from  the  difference  in  the  positions  of  the  two  extrema  {ymaz->  ymin^ 

=  ytnax  -  ymin-  (3.4.2  o) 

The  maximum  value  of  the  streamwise  vorticity  intensity  is  Ein  estimate  of  the 
strength  (7)  of  the  Rankine  vortex,  and  the  tangential  velocity  at  the  edge  of  the 
vortex  is  given  by 

vt  =  7*r+.  (3.4.2  S) 


The  location  of  the  maiximum  of  the  streamwise  intensity  {ytnax)  ^  estimate 
of  the  mean  position  of  the  center  of  the  vortex  core.  From  the  circulation  of  the 
near-wall  Rankine  vortex 


2ir 

j  VeVedB, 

0 


(3.4.2  c) 
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a  vortex  circulation  Reynolds  number  can  be  defined  by 


r 

jRcf  =  —  =  27rT;Jrg .  (3.4.2  d) 

These  parameters  of  the  near- wall  streamwise  vortices  are  listed  in  Table  3.3  for 
the  two  transversely  curved  flows  as  well  as  the  plane  channel  flow  of  Kim  et  al. 
(1987).  It  is  noteworthy  that  neither  the  core  radius  nor  the  position  of  the  center 
of  these  vortices  changes  appreciably  with  curvature,  even  though  the  strength  of 
the  Rankine  vortex  is  greatly  reduced. 

There  axe  other  important  curvatme  effects  on  the  vorticity  intensities.  Near 
the  wall  (y*"  <  7),  the  normal  vorticity  intensity  is  virtually  unaffected  by 

the  cmrvature,  while  the  other  vorticity  components  decrease.  In  the  outer  part  of 
the  flow  as  shown  in  Figure  3.17. 

The  inclination  of  the  projection  of  the  vorticity  vector  in  (r,  2:)-planes  is  given 

by 


0^2  =  ^ 


(3.4.3) 


The  probabiUty  density  functions  (p.d.f.s)  of  0rr,  weighted  by  the  magnitude  of 
the  the  projected  vorticity  vector  (Moin  &  Kim  [1985]),  (w*  +^l)/<  <^r  +<^z  >> 
axe  shown  in  Figures  3.18  o  and  3.18  b.  In  the  expression  above,  <>  indicates 
the  rnesm  of  the  quantity  inside  the  brackets  taken  on  the  corresponding  {z,6) 
cylindric2J  surface.  The  following  discussion  refers  to  the  weighted  p.d.f.s,  even 
though  all  the  features  described  below  are  also  evident  in  the  unweighted  p.d.f.s. 
The  advantage  of  the  weighting  is  that  it  enhances  the  contributions  of  the  strong 
vorticity  fluctuations. 

Throughout  the  discussion  of  the  orientation  of  the  vorticity  vector,  the  reader 
should  bear  in  mind  that  the  plane  channel  flow  fields  to  which  these  results  axe 
compared  were  from  the  large  eddy  simulation  of  Moin  &  Kim  (1985),  which  was 
performed  for  a  Reynolds  number  of  13800  (based  on  the  centerline  velocity  and 
^).  At  the  closest  point  to  the  wall  {y*  w  4)  invest:^^  ted  by  Moin  &  Kim  (1985), 
the  distribution  of  Qrz  is  centered  around  0®  and  ±180®  since  the  normal  vorticity 
must  go  to  zero  at  the  wall  due  to  no-slip.  In  the  transversely  curved  flows,  this 
effect  is  not  discernable  in  the  histograms  of  0rz  for  y"*"  >  2.  At  y*"  %  4,  both 


5 


3- . 

coordinates  fo  ^  ^  J  ^  intensity. 
a1  intensity - iNoriu<« 


Figure  3.18  a  Distribution  of  the  inclination  of  the  angle,  ©r*,  of  the  projection  of  the  vorticity  vectors  in 
(r,  2:)-planes  for  6/a  =  5.0;  data  weighted  with  the  magnitude  of  the  projected  vorticity.  (a)  y'*'  =  0.66;  (6) 
y+  =  4.10;  (c)  y+  =  8.02;  (d)  y+  =  16.31;  (e)  y+  =  25.38;  (f)  y+  =  38.76;  (o)  y+  =  51.85;  (h)  y+  =  69.66;  (i) 
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planes  for  Sfa  —  11.0;  data  weighted  with  the  magnitude  of  the  projected  vorticity. 
4.46;  (c)  =  8.74;  (d)  =  17.77;  (e)  y+  =  27.66;  (/)  y+  =  42.24;  (y)  y+  =  56.52; 

104.90. 


curved  flows  have  distributions  of  ©rz  that  peak  around  ±90®.  This  distribution  of 
0rz  persists  up  to  y'*’  ss  25  in  the  Sfa  =  5  flow  and  up  to  y"*"  w  40  in  the  S/a  =  11 
flow.  Father  from  the  wall,  the  the  distributions  gradually  broaden,  and  by  y'*'  «  52 
in  the  Sja  =  b  flow  and  y'*'  »  75  in  the  S/a  =  11  the  peak  shifts  to  —135®  and  45® 
as  was  observed  by  Moin  &  Kim  (1985).  It  is  only  beyond  these  points,  which  is 
most  of  the  domain,  that  the  concentration  of  the  distributions  of  Qrz  shift  to  the 
—135®  and  45®  orientations. 

In  Moin  &  Kim  (1985)  the  main  mechanism  of  vorticity  stretching  in  the  plane 
channel  was  attributed  to  stretching  by  the  mean  shear,  which  has  its  principal  axis 
at  45®  with  the  direction  of  the  mean  flow.  They  also  pointed  out  that,  according 
to  Deissler  (1969),  when  the  normal  and  streamwise  vorticity  intensities  are  equal, 
the  direction  of  the  majcimum  vorticity  in  a  shear  flow  is  at  45®  to  the  mean  flow. 
In  the  transversely  curved  flows  there  is  an  increa-^ingly  thicker  laj'er  around  the 
cylinder,  in  which  the  vorticity  has  the  ±90®  orientation.  This  suggests  that,  as 
the  curvature  increases,  vortical  structures  inclined  at  45®  to  the  mean  flow  are 
probably  more  difficult  to  observe. 

3.5  Instantaneous  Turbulent  Flow  Structures 

Contours  of  the  streamwise  velocity  fluctuations  (vz)  on  two  (r,  $)  planes  (normal 
to  the  mean  flow)  eire  shown  in  Figures  3.19  (plane  Iz)  and  3.20  (plane  IIz).  In  the 
smaller  curvature  case  {S/a  =  5,  see  Figures  3.19  a  and  3.20  o),  the  perimeter  of 
the  cylinder  is  about  270  wail  units  and  in  the  snap  shots  shown  only  three  to  four 
low-speed  streaks  are  observed  around  the  cylinder.  In  the  strongly  curved  case 
{S/a  =  11,  see  Figures  3.19  b  and  3.20  6),  the  perimeter  of  about  140  wall  units 
apparently  can  only  support  two  low-speed  streaks.  A  better  measure  of  the  mean 
spacing  between  the  low-speed  streaks  is  given  by  the  velocity  correlations  in  the 
azimuthal  direction  at  y"*"  w  12,  shown  in  Figure  3.21.  The  mean  streak  spacing  is 
about  100  wall  units  in  both  flows  at  this  radiail  location.  Note  however  that  the 
mean  streak  spacing  has  a  nearly  linear  dependence  on  the  distance  from  the  wall 
(y"*")  with  a  slope  that  increases  with  curvature  (Figure  3.22). 

Figures  3.23  and  3.24  show  contours  of  the  radial  (wall-normal)  velocity  fluc¬ 
tuations  in  the  same  planes  as  Figures  3.19  and  3.20  respectively.  The  constraint 
imposed  by  the  outer  boundary  conditions  (Equation  2.1.5  b)  on  the  radial  motion 
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Figure  3.19  Contours  of  axial  velocity  fluctuations  (uz),  normalized  by  Ur,  on  a 
plane  {Iz)  normal  to  the  mean  velocity:  (a)  6 /a  =  5  with  contour  levels  from  —5.2 
Ur  to  3.55  Ur;  (6)  8/a  =  11  with  contour  levels  from  —4.95  Ur  to  4.05  Ur-  The 
contour  increments  of  0.25  Ur-  The  solid  contours  denote  the  low-speed  fluctuations 
{vz  <  0)  and  the  dotted  contours  denote  the  high-speed  fluctuations  {vz  >  0). 


56 


Figure  3.20  Contours  of  axial  velocity  fluctuations  (u;),  normalized  by  Ur,  on  a 
plane  {Hz)  normal  to  the  mean  velocity:  (a)  6/a  =  5  with  contour  levels  from  —5.2 
Ur  to  3.3  Uri  {h)  6 ! a  =  W  with  contour  levels  from  —4.2  Ur  to  3.3  Uj-  The  contour 
increment  is  0.25  Ut-  The  solid  contours  denote  the  low-speed  fluctuations  {  vz  <  0) 
and  the  dotted  contours  denote  the  high-speed  fluctuations  {vz  >  0). 


Figure  3.21  Azimuthal  two-point  correlations  of  the  velocity  fluctuations  for  (a) 
6/a  =  5.0  and  for  (6)  b/a  =  11:  - - Rv^Vr^ - Rvgve- 
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Figure  3.22  The  mean  spanwise  streak  spacing  A"*"  estimated  from  the  two-point 
correlations  of  -t-  plane  channel  (Kim  tt  al.  [1987]);  cylinders:  o,  for  6/a  =  5 
and  A  ,  for  S/a  =  11. 

of  structures  in  the  outer  part  of  the  layer  is  evident  particularly  in  the  6 /a  =  5  flow 
(Figures  3.23  a  and  3.24  a).  From  the  contour  plots  of  Vz  and  Vr,  regions  of  sweep 
2ind  ejection  of  fluid  to  and  away  from  the  cylinder  surface  can  be  easily  identified. 

In  simulations  of  homogeneous  turbulence  with  a  small  shear  rate,  Rogers  & 
Moin  (1987)  showed  that,  even  though  the  contours  of  the  streamwsie  velocity 
fluctuations  tend  to  be  elongated  in  the  streamwise  direction,  no  strezdcy  structures 
were  observed.  In  studies  of  homogeneous  turbulence  at  a  high  shear  rate,  Lee  et 
al.  (1990)  observed  long  streaks  in  the  absence  of  a  no-slip  wall,  and  concluded 
that  it  is  the  magnitude  of  the  mean  shear  rate  that  determines  whether  streaks 
are  observed.  These  two  studies  suggest  that  the  streamwise  extent  of  the  streaks 
is  controlled  by  the  mean  shear.  A  measure  of  the  strength  of  the  mean  shear  is 
the  shear  rate  parameter  (Lee  et  al.  [1990]) 
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Figure  3.23  Contours  of  radial  (wall- normal)  velocity  fluctuations  (vr),  normalized 
by  Ur,  on  a  plane  {Iz)  normal  to  the  mean  velocity:  (a)  Sja  =  b  with  contour  levels 
from  —1.8  Ur  to  1.3  Ur;  (b)  6/a  =  11  with  contour  levels  from  —0.95  ur  to  1.05 
Ur.  The  contour  increment  is  0.25  Ur-  The  solid  contours  denote  the  negative 
radial  velocity  (vr  <  0)  and  the  dotted  contours  denote  the  positive  radial  velocity 
(ur  >  0). 
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Figure  3.24  Contours  of  radial  (wall-normal)  velocity  fluctuations  (ur),  normalized 
by  Ur,  on  a  plane  {IIz)  normal  to  the  mean  velocity:  (a)  Sfa  =  5  with  contour 
levels  from  —1.05  Ur  to  1.3  Ur;  (b)  6 /a  =  11  with  contour  levels  from  —0.95  Ur  to 
1.29  Ur-  The  contour  increment  is  0.25  Ur-  The  solid  contours  denote  the  negative 
radial  velocity  (vr  <  0)  and  the  dotted  contours  denote  the  positive  radial  velocity 
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from  2.2  Ur  to  4.85  Ur-  The  axial  length  of  the  domain  displayed  is  and  the  contour  increment  is  0.25  Ut-. 
The  solid  contours  denote  the  low-speed  fluctuations  {vz  <  0)  and  the  dotted  contours  denote  the  high-speed 
fluctuations  (vz  >  0). 


where  and  are  the  ajciaJ  and  azimuthal  integral  scales  of  the  quantity  of 
interest.  For  example,  when  based  on  and  are  given  by 


h,l‘l 

^v^(y)=  j 

0 


LvKv)  =  («  +  y)  J 

0 


(3.5.3) 


respectively.  L*  based  on  the  streamwise  velocity  fluctuations,  which  characterize 
the  low  speed  streaks,  was  computed  by  Lee  et  al.  (1990)  for  the  plane  channel; 
their  results  are  compared  to  the  present  simulations  in  Figure  3.27.  Close  to  the 
wall  {y'^  <  10),  L*  is  signincantly  increased  with  increzising  curvature.  Since  the 
spanwise  length  scale  of  the  low-speed  streaks  in  wall  units  is  not  greatly  affected 
by  curvature,  their  streamwise  length  sccJe  must  increase  with  6/ a.  Away  from  the 
wall  (j/+  >  20),  L*  for  the  S/a  =  h  cylinder  has  essentially  the  same  value  as  in 
the  plane  channel.  In  the  larger  curvature  case  {S/a  =  11)  L*  away  from  the  wall 
(y"*"  >  20)  is  twice  as  large  as  its  counterparts  for  the  S/a  =  5  cylinder  and  the 
plane  channel. 

In  addition  to  the  low  speed  streaiks,  near-wall  flows  are  also  characterized  by 
neair-wall  vortical  structures  and  internal  shear  layers.  The  near- wall  streamwise 
vortices  and  internal  shear  layers  are  characterized  by  the  streamwise  aind  the  span- 
wise  vorticity  fluctuations,  and  w#,  respectively.  The  length  scale  ratios  defined 
in  Equation  3.5.2,  computed  for  both  Wj  aind  a;#,  are  shown  in  Figure  3.28.  For 
J/'*'  >  10,  the  length  scale  ratios  for  the  two  vorticity  components  are  close  to  their 
isotropic  value  of  unity.  Near  the  wall  (y'*'  <  10),  where  the  mean  shear  is  lairgest 
(Figure  3.7),  the  length  scale  ratios  of  ug  and  u^z  increase  with  increasing  curva¬ 
ture.  In  particular,  the  higher  neair-wall  L*  based  on  ujz  suggests  that  the  near-wall 
streamwise  vortices  become  longer  as  the  curvature  increases.  Note  also  that  both 
vorticity  length  scale  ratios  are  only  affected  by  curvature  near  the  wall  (y'*'  <  10). 

The  near- wall  shear  layers  that  develop  at  the  interface  of  the  low  and  high  speed 
flow  regions  are  common  in  the  two  transversely  curved  flows  as  shown  in  Figures 
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Figure  3.27  The  length  scale  parameter  L*  based  on  Vz\ - plane  channel  (Lee 

ei  al.  [1990]);  cylinders  with - S/a  =  5  and - 6/a  =  11. 


q 


Figure  3.28  The  length  scede  parzuneter:  L*  based  on  for  cylinders  with 

- S/a  =  5  2uid - S/a  =  11;  L*  based  on  ujz  for  cylinders  with - S/a  =  5 

emd .  S/a  =  11. 
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Figure  3.29  Contours  of  axial  (streamwise)  vorticity  fluctuations  normalized 
by  Ur  and  u,  on  a  plane  {Iz)  normal  to  the  mean  velocity;  (a)  5/a  =  5  with  contour 
levels  from  -0.34  u\lv  to  0.5  u\lv\  (6)  5/a  =11  with  contour  levels  from  -0.3 
u\lv  to  0.38  u\fv.  The  contour  increment  is  0.04  u^/u.  The  solid  contours  denote 
the  negative  vorticity  {uz  <  0)  and  the  dotted  contours  denote  the  positive  vorticity 
[Uz  >  0). 
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(b) 


Figure  3.30  Contours  of  axial  (stresimwise)  vorticity  fluctuations  normalized 
by  Ur  and  u,  on  a  plane  (IIz)  normal  to  the  mean  velocity;  (a)  6 /a  =  5  with 
contour  levels  from  —0.34  u^/i/  to  0.5  u\lv\  (6)  S/a  =  11  with  contour  levels  from 
—0.26  ui/i/  to  0.22  u\lv.  The  contour  increment  is  0.04  The  solid  contours 

denote  the  negative  vorticity  (uz  <  0)  and  the  dotted  contours  denote  the  positive 
vorticity  {u>z  >0). 
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Figure  3.33  Contours  of  azimuthal  vorticity  fluctuations  (u;,),  normalized  by  Uj 
and  u,  on  a  plane  (Iz)  normal  to  the  mean  velocity;  (a)  6/a  =  5  with  contour 
levels  from  -0.58  u\/u  to  0.54  (b)  6/a  =  11  with  contour  levels  from  -0.66 

/v  to  0.5  ul/u.  The  contour  incremen*^^  is  0.04  u\/v.  The  solid  contours  denote 
the  negative  vorticity  (w,  <  0)  (aligned  with  the  mean  vorticity)  and  the  dotted 
contours  denote  the  positive  vorticity  {ug  >  0). 
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Figure  3.34  Contours  of  azimuthal  vorticity  fluctuations  (u;#),  normalized  by  Ur 
and  I/,  on  a  plane  (IIz)  normal  to  the  mean  velocity:  (a)  S/a  =  5  with  contour 
levels  from  —0.38  u^/u  to  0.5  u\jv\  (b)  Sja  =  11  with  contour  levels  from  —0.9 
u\lv  to  0.46  u\lv.  The  contour  increment  is  0.04  u\lv.  The  solid  contours  denote 
the  negative  vorticity  (w#  <  0)  (aligned  with  the  mean  vorticity)  and  the  dotted 
contours  denote  the  positive  vorticity  (u;#  >  0). 


levels  from  -0.58  u^i/  to  0.32  (6)  <5/a  =  11  with  contour  levels  from  -0.73  u\lv  to  0.62  u\lv.  The 

contour  increment  is  0.05  u\lv.  The  solid  contours  denote  the  negative  vorticity  <  0)  (aligned  with  the 
mean  vorticity)  and  the  dotted  contours  denote  the  positive  vorticity  (u;«  >  0). 


3.31  (plane  Iz)  and  3.32  (plane  IIz)  and  have  features  similar  to  those  of  the  plane 
channel.  These  two  figures  suggest  that  the  higher  curvature  {6 /a  =  11)  flow  is 
more  quiescent  in  the  outer  region  in  agreement  with  the  lower  velocity  correlation 
coeflBcient  (Figure  3.6).  In  Figures  3.33  and  3.34,  cross  sections  of  the  internal  shear 
layers  by  (r,  0)  planes  show  that  the  shear  layers  have  large  azimuthal  (spanwise) 
length  scales  relative  to  the  cylinder  radius.  As  the  curvature  increases  the  ratio 
of  the  spanwise  length  scale  of  the  shear  layers  to  the  cylinder  radius  increases. 
These  contour  plots  also  suggest  that  the  shear  layers  have  a  milder  transverse 
curvature  than  the  cylinder.  Cuts  of  these  shear  layers  by  (z,  9)  surfaces  parallel 
to  the  cylinder  (Figure  3.35)  at  y'*'  «  15  show  a  characteristic  arrow  shape  that  is 
more  noticeable  in  the  higher  curvature  (6/ a  =  11)  flow.  Such  arrow  shapes  are 
geometrically  consistent  with  nearly  planar  inclined  shear  layers  intersected  by  a 
cylindrical  surface. 

Kim  (1989)  demonstrated  that  the  near-wall  streamwise  vortices  contribute  sig¬ 
nificantly  to  the  source  term  of  the  pret-  are  Poisson  equation.  As  the  curvature 
increases,  there  is  an  additional  important  source  of  pressure  fluctuations  (see  Sec¬ 
tion  4.1,  Figure  4.3)  associated  with  the  strong  near-wall  radial  (wall-normal)  vor- 
ticity  fluctuations  (Figures  3.18  a  and  3,18  b).  Contours  of  the  radial  vorticity 
fluctuations  in  the  same  {r,0)  planes  are  shown  in  Figures  3.36  and  3.37.  As  in 
the  plane  channel,  near  the  wall  the  regions  where  the  radial  vorticity  fluctuations 
change  sign  are  well  correlated  with  the  locations  of  the  low-speed  streaks  but  not 
necessarily  with  the  streamwise  vortices.  Contours  of  the  pressure  in  the  same 
(r,  0)  planes  (Figures  3.38  and  3.39)  show  low  pressure  regions  well  correlated  with 
the  locations  of  the  near-wall  streamwise  vortices.  Apparently  the  normal  vorticity 
fluctuations  become  a  more  prominent  feature  simply  because,  as  the  curvature 
increases,  more  of  the  turbulent  kinetic  energy  is  in  the  streamwise  velocity  fluc¬ 
tuations  (Figure  3.4).  Contours  of  the  normal  vorticity  ^lre  shown  in  Figure  3.40 
on  {z,0)  surfaces  parallel  to  the  cylinder  at  «  15.  In  the  6/a  =  11  flow  (Figure 
3.40  b),  the  strong  posicive  and  negative  contours  of  u)r  eire  well  correlated  with  the 
intersections  of  the  cylindrical  surface  at  y*  «  15  with  the  wall  shear  layers  (Figure 
3.35  b). 
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Figure  3.36  Contours  of  radial  (wall-normal)  vorticity  fluctuations  (wr),  normalized 
by  Ur  and  v,  on  a  plane  (/r)  normal  to  the  mean  velocity:  (a)  6/a  =  5  with  contour 
levels  from  -0.38  ul/v  to  0.46  (b)  8/a  =  11  with  contour  levels  from  -0.3 

Ur/t'  to  0.38  u\lv.  The  contour  increment  is  0.04  u\lv.  The  solid  contours  denote 
the  negative  vorticity  {ujr  <  0)  and  the  dotted  contours  denote  the  positive  vorticity 
(wr  >  0). 
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F/gure  3.37  Contours  of  radial  (wedl-normal)  vorticity  fluctuations  normalized 
by  Ur  and  t/,  on  a  plwe  (IIz)  normal  to  the  meeui  velocity:  (a)  S /a  =  5  with  contour 
levels  from  —0.46  u\lv  to  0.5  u\lv\  (6)  S/a  =  11  with  contour  levels  from  —0.26 
u\/u  to  0.22  u\/u.  The  contom-  increment  is  0.04  u\/v.  The  solid  contours  denote 
the  negative  vorticity  (a>r  <  0)  and  the  dotted  contours  denote  the  positive  vorticity 
(wr  >  0). 
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Figure  3.38  Contours  of  the  pressure  fluctuations  (p),  normalized  by  p  and  Ur, 
on  a  plane  (Iz)  normal  to  the  mean  velocity:  (a)  6/a  =  5  with  contour  levels  from 
—4.1  pu%  to  2.1  pu’;  (6)  6/a  =  11  with  contour  levels  from  —2.5  pu%  to  1.5  pu%. 
The  contour  increment  is  0.2  pu^.  The  solid  contours  denote  the  negative  pressure 
(p  <  0)  and  the  dotted  contours  denote  the  positive  pressure  (p  >  0^ 
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Figure  3.39  Contours  of  the  pressure  fluctuations  (p),  normalized  by  p  and  Ur?  on 
a  plane  {IIz)  normal  to  the  mean  velocity:  (a)  6 /a  =  5  with  contour  levels  from 
—2.5  pu^  to  1.5  (b)  S/a  =  11  with  contour  levels  from  —1.7  pu^  to  0.5  pu^. 

The  contour  increment  is  0.2  pu^.  The  solid  contours  denote  the  negative  pressure 
(p  <  0)  and  the  dotted  contours  denote  the  positive  pressure  (p  >  0). 
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CHAPTER  4 


Pressure  Fluctuations 


Wall  pressure  fluctuations  affect  the  acoustic  characteristics  (both  radiation  and 
detection)  of  underwater  vehicles.  In  addition,  in  transversely  curved  turbulent 
flows  over  sonar  devices,  it  is  important  to  identify  the  pressure  signals  generated 
by  the  turbulent  flow  field.  The  aim  of  this  chapter  is  to  describe  the  effects  of 
transverse  curvature  on  the  wall  pressure  fluctuations  when  both  the  outer  and 
inner  parts  of  the  flow  are  affected  by  the  curvature  of  the  wall.  For  this  flow 
regime  (large  S/a  and  small  a"*'),  there  is  only  limited  experimental  data  available 
on  the  statistical  properties  of  the  wall  pressure  fluctuations. 

For  compeurison,  the  spatial  and  temporal  spectra  of  the  plane  chainnel  wall 
pressure  fluctuations  were  also  calculated  from  a  database  generated  by  Kim  et 
al.  (1987)  and  Choi  &:  Moin  (1990).  As  noted  earlier,  the  Reynolds  number  of  the 
plane  channel  flow  and  those  of  the  present  flows  are  similar,  and  thus  differences 
in  the  statistical  characteristics  can  be  ascribed  to  the  transverse  curvature  alone. 

The  wall  pressure  database  on  which  the  results  of  this  chapter  are  based  was 
generated  after  the  flow  had  reached  statisticed  steady  state  (Section  3.2).  Addi¬ 
tional  integration  of  the  governing  equations  for  wall  pressure  statistics  were  carried 
out  for  a  period  of  12.8  S/ut  time  units  in  the  5/a  =  5  flow  and  for  8.6  time  units 
in  the  5/a  =  11  flow. 

The  effects  of  transverse  curvature  on  the  sources  of  pressure  fluctuations  are 
discussed  in  Section  4.1.  The  computational  method  used  to  obtain  the  space-time 
characteristics  of  the  wall  pressure  data  is  described  in  Section  4.2.  Section  4.3 
contains  the  spatial  spectra  of  the  wall  pressure  fluctuations  and  in  Section  4.4  the 
two-point  correlations  are  presented.  The  fractional  contributions  of  various  flow 
regions  to  the  wall  pressure  are  studied  in  Section  4.5.  The  temporal  spectra  are 
presented  in  Section  4.6.  Two-dimensional  spectra  and  correlations  are  discussed  in 
Section  4.7.  Finally,  Section  4.8  treats  the  convection  velocity  of  the  wall  pressure 
fluctuations. 
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Figure  4.1  Root-mean-squaxe  pressure  normalized  by  Ur  in  global  coordinates: 
- plane  channel  (Kim  [1989]);  cylinders  with - 6/a  =  5  and - 8/a  =  11. 

4.1  Pressure  Source  Terms 

The  Poisson  equation  for  the  pressure  field  is  obtained  from  the  divergence  of 
the  Navier-Stokes  equations  (Equation  2.1.1), 


^dVzdvr 
dr  dz 


dvz  dvr 
dr  dz 


2  ( dvr  \ 

r  dr  \dd 

2  dvz  dvg  /  dvz  \  ^ 
r  d0  dz  \dz  ) 


(4.1.1  a) 


Boundary  conditions  for  the  pressure  are  obtained  by  evaluating  the  radial  momen¬ 
tum  equation  at  both  edges  of  the  domain, 


dr 


1  d^r 

dp 

1  2  dvg  vg^' 

Rer  dr^ 

1 

r=a 

dr 

r=a-|-l 

Rer  r^  do  ^  r 

r=a-t-l 


(4.1.1  b) 
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Figure  4.2  Profiles  of  the  root-mean-square  value  of  the  pressure  source  terms 

normalized  by  v  and  Ut'-  (a)  linear  contribution:  cylinders  with  -  6/a  =  ^ 

and - 6/ a  =  11;  (6)  nonlinezu"  contribution:  cylinders  with - 6/a  =  5  and 

- 6 /a  =  11;  total  source  terms:  cylinders  with - 6/a  =  5  and - 6/a  =  11. 
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The  first  term  on  the  right  hand  side  of  Equation  4.1.1  a  represents  the  pressure 
source  that  results  from  the  linear  interaction  of  the  mean  shear  with  the  turbulence. 
The  next  six  terms  are  the  pressure  source  terms  that  result  from  the  nonlinear 
interactions  within  the  turbulent  flow  field. 

The  pressure  intensity  normalized  by  the  mean  shear,  shown  in  Figure  4.1,  de¬ 
creases  across  the  layer  with  increasing  curvature.  In  particular,  the  near-wall 
maximum  of  the  pressure  intensity  is  reduced  by  increasing  the  curvature.  This 
suggests  that  the  structures  that  are  responsible  for  the  maximum  in  the  pressure 
intensity  are  weakened. 

Contrary  to  what  was  generally  accepted,  Kim  (1989)  found  that  in  turbulent 
channel  flow  the  turbulence  field  interactions  constitute  the  strongest  sources  of 
pressure  fluctuations.  In  fact,  the  magnitude  of  the  mean-square  of  the  linear 
source  term  was  about  five  times  lower  than  the  magnitude  of  the  mean-square 
of  the  nonlinear  source  term.  In  the  transversely  curved  turbulent  flows  reported 
here,  the  ratio  of  the  magnitude  of  the  mean-square  of  the  nonlinear  to  the  linear 
pressure  source  terms  is  large  (Figure  4.2),  even  though,  as  the  curvature  increases, 
both  sources  of  pressure  fluctuations  (linear  and  nonlinear)  decrease. 

In  the  plane  channel  the  maximum  of  the  mean  square  of  the  nonlinear  terms 
occurs  at  about  «  20,  which  is  the  same  as  the  mean  position  of  the  near-wall 
vortices  as  well  as  the  position  of  the  maximum  of  the  pressure  intensity  (Kim 
[1989]).  In  the  transversely  curved  flows,  the  maximum  of  the  pressure  intensity  is 
weakened  as  the  curvature  increases  and  at  about  the  same  position  as  in  the  plane 
channel  flow  (j/’*’  20,  Figure  4.1).  As  the  curvature  increases,  the  average  strength 

of  the  near- wall  vortices  is  weakened  while  the  average  position  of  the  vortex  cores 
remains  the  same  (y"*"  w  20)  (see  Section  3.4). 

In  order  to  associate  the  pressure  intensity  with  the  sources  of  pressure  fluctu¬ 
ations  in  the  flow,  the  various  contributions  to  the  nonlinear  source  terms  of  the 
pressure  are  examined  (Figure  4.3).  As  in  the  plane  channel  the  most  important 
nonlinear  source  term  is 


which  contains  significant  contributions  from  the  neeir-wall  streamwise  vortices. 
Note  also  that  in  both  flows  this  source  term  develops  near-W2ill  (y"*"  sa  5)  local 
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Figure  4.3  Profiles  of  the  root-mean-square  value  of  the  nonlinear  contributions  to 

the  source  terms  normalized  by  u  and  Uj:  (a)  b/a  =  5,  (6)  b/a  =  11; - {dvr/drY, 

.  {2lr){dvrlde-v,){dv,ldT), - (l/^2)(d^;,/a^?-t-Ur)^  - {dv,ldz)\  — - 

2{dv,ldr){dvrldz), - {2lr){dv,lde){dv,ldz). 
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Lt{'^tI8) 

9.36 

12.9 

8.6 

Ltiur/ 8) 

1.20 

0.806 

0.659 

Ats(ur/8) 

3.75  X  10"^ 

2.52  X  10“^ 

2.06  X  10 

Ats{ullu) 

0.675 

0.539 

0.494 

A(jj{8/ut) 

5.2 

7.8 

9.5 

Ns 

2560 

5120 

4160 

Table  4.1  Temporal  resolution  parameters 


extrema.  Near  the  wall  vorticity  of  opposite  sign  to  the  primary  vortices  may  be 
responsible  for  this  behavior. 

Figure  4.3  also  shows  that,  as  the  curvature  increases,  another  nonlinear  term. 


2  dvz  dvg 


(4,1.3) 


plays  an  increasingly  important  role  in  the  pressure  fluctuations.  Note  that  the 
factors  in  this  product  are  the  two  velocity  gradients  that  define  the  radial  compo¬ 
nent  of  the  vorticity  (Equation  3.4.1).  As  the  curvature  increases,  strong  normal 
vorticity  fluctuations  (with  respect  to  the  local  streamwise  vorticity  fluctuations) 
become  increasingly  more  common  farther  away  from  the  wall  (when  measiured  in 
wall  units). 


4.2  Computational  Considerations 

The  wall  pressure  was  sampled  at  intervals  At®  for  a  total  period  Lj’,  resulting 
in  Ng  time  samples.  In  the  temporal  analysis  a  window  of  length  Lt  is  used.  These 
pMameters  are  given  in  Table  4.1  for  the  three  simulations  under  consideration. 
For  the  sake  of  simplicity  only  the  one-dimensional  spectra  are  discussed.  The 
computation  of  the  two-dimensional  spectra  is  identical. 
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4.2.1  Spatial  Spectra 

The  spatial  spectra  of  the  wall  pressure  fluctuations  are  obtained  by  averaging 
over  all  Ns  samples.  Thus  the  spectra  in  the  axial  and  azimuthal  directions  are 
given  by 


Epp{kz) 


N, 

E 

i=i 


■  Ns/2 

Pw{kz,Tn6/a,tj)pl,{kz,Tn6/a,tj)  , 

m=-Ns/2 


Epp{ks)  — 


^ksN, 


N. 

E 

i=i 


X/2 

^  Pwi2irk/Lz,ks,tj)pl,{2irklLz,ks,tj)  , 
Jfc=0 


(4.2.1) 


respectively.  The  Fourier  transform  of  pw{z,0,tj)  is  pw{kz,  ks,tj)  and  pXo{kz,  ks,ij) 
is  its  complex  conjugate.  Note  that  the  mean-square  wall  pressure  is  given  by 


N^/2  Ns/2 

p^yj  =  Akz  ^  Epp(27rk/Lz)  =  Aks  ^  Epp{mS/a).  (4.2.2) 

jb=0  m=0 


4.2.2  Temporal  Spectra 

In  the  cedculation  of  the  temporal  characteristics  of  the  wall  pressure  the  same 
techniques  were  used  as  in  Choi  &  Moin  (1990).  The  total  length  of  the  temporal 
domain  for  each  of  the  flows  (Tj'),  is  divided  into  m  overlapping  intervals  (with 
50%  overlap)  of  equal  length  (Tt).  The  length  of  each  interval  is  thus  given  by 


Lt  = 


2Lj’ 
m  4- 1 


(4.2.3) 


The  choice  of  m  =  31  and  m  =  25  was  made  for  the  S/a  =  5  and  S/a  =  11  flows 
respectively,  resulting  in  320  samples  in  each  interval.  Since  the  wall  pressure  is 
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not  temporally  periodic  in  any  of  the  intervals,  the  data  is  tapered  at  the  edges  of 
each  interval  by  a  window  function  fwii),  according  to 


Pw{t)  =  fw{t)pw{t),  0<t<Lt,  (4.2.4) 

where  pw{t)  denotes  the  wall  pressure  at  a  given  time  instant.  For  the  sahe  of 
simplicity,  in  this  section,  the  spatial  dependence  of  the  wall  pressure  fluctuations 
is  omitted.  In  the  anzdysis  of  the  wall  pressure  for  the  plane  channel  Choi  &  Moin 
(1990)  experimented  with  several  window  functions  and  reported  no  signiflcant 
difference  in  the  resulting  temporal  spectra.  In  the  present  calculations  the  Hanning 
window  function 

...  1  1  ( 2nt\ 

Mt)  =  2  “  2  \m)  ’ 

is  used.  The  properties  of  fw(i)  as  a  window  function  are  discussed  by  Harris 
(1978). 

The  temporal  spectrum  is  estimated  from 

<f>pp(uj)  =  (4-2.6) 


where  p^,  is  the  Fourier  transform  of  pw  in  the  interval.  The  discrete  frequencies, 
u)j,  axe  given  by 

j  =  0,1, 2,...,  M/2,  (4.2.7) 

where  M  =  2Ns/{m  -f  1)  is  the  number  of  samples  in  the  interval.  In  general, 
because  of  the  windowing,  the  variance  of  pu,  will  not  be  equal  to  the  variance  of 
Pv).  To  preserve  the  root-mean-squaxe  (r.m.s.)  value  of  pw  the  spectrum  <f)pp  is 
rescaled  so  that 

M/2 

Pw  =  <f>pp{^j)-  (4.2.8) 

j=-M/2 
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Figure  4.4  Signal  of  the  wall  pressure  fluctuations  normalized  by  : 
channel;  cylinders  with  (6)  6/a  =  5;  (c)  S/a  =  11. 


(a)  plane 
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Sfa 

0 

5 

11 

¥'wl{pu\) 

1.54 

0.84 

0.58 

S{pw) 

-0.03 

-0.70 

-0.78 

npw) 

5.16 

7.16 

6.18 

Table  4.2 

Wall  pressure  parameters. 

4.3  Spatial  Spectra 

Figure  4.4  shows  a  sample  of  the  time  traces  of  the  wall  pressure  fluctuations 
normalized  by  their  r.m.s.  values.  Note  that,  as  the  curvature  increases,  the  high 
frequency  content  of  the  fluctuations  is  weakened.  Apparently  the  different  appear¬ 
ances  of  the  wall  pressure  fluctuations  in  the  two  transversely  curved  flows,  does 
not  noticeably  affect  the  flatness  and  skewness  values  (Table  4.2). 

The  strong  dependence  of  the  wall  pressure  intensity  on  the  Reynolds  number, 
documented  by  Choi  &  Moin  (1990),  is  not  a  factor  in  the  comparison  of  Table  4.1 
because  the  Reynolds  numbers  of  the  three  simulations  are  similar.  In  measure¬ 
ments  of  the  wall  pressure  intensity  in  a  boundary  layer  Willmarth  &  Yang  (1970) 
concluded  that  the  transverse  curvature  did  not  have  a  significant  effect.  However, 
this  result  is  for  a  mild  transverse  cxirvature  (^/a  ~  2,  a+  ~  4500),  in  which  only 
the  outer  part  of  the  boundary  layer  was  affected  by  the  curvature.  In  addition,  any 
small  curvature  dependence  that  was  present  in  their  study  was  probably  overshad¬ 
owed  by  the  strong  Reynolds  number  dependence  of  the  wall  pressure  intensity.  In 
subsequent  wzdl  pressure  measurements  on  a  cylinder  (Willmarth  et  al.  [1976]),  for 
which  6/a  ~  4,  the  value  of  the  wall  r.m.s.  pressiu'e  was  not  reported. 

The  one- dimensional  waJl  pressure  spectra  as  a  function  of  the  azimuthal  (span- 
wise)  Euad  axial  (streamwise)  wave-numbers  are  shown  in  Figure  4.5,  along  with 
their  planar  counterparts  (Choi  &  Moin  [1990]).  Note  that  the  spzinwise  spectra 
(Figure  4.5  o)  is  not  affected  by  the  transverse  cinrvature  in  the  high  wavenumber 
range.  In  the  low  wavenumber  range  the  energy  of  the  spanwise  spectrum  decreases 
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Figure  4.6  Axial  (streamwise)  spectra  of  the  wall  pressure  fluctuations  nonnedized 

by  6j  and  Ur". - plane  channel  (Choi  &  Moin  [1990])  {6 fa  =  0,  =  6)\  cylinders 

with - 8 la  =  5  and - 8 fa  =  11. 

with  increasing  curvature.  This  suggests  that  the  curvature  mostly  affects  struc¬ 
tures  with  large  spanwise  length  scales  (small  8kt)  than  those  with  smedl  spanwise 
length  scales  (large  8k$).  The  streamwise  spectral  density  of  the  wall  pressure 
fluctuations  decreases  significantly  with  curvature  for  all  scales  (Figure  4.5  6).  As 
expected,  the  three  flows  display  a  negligible  wave-number  range  with  the  —1  slope 
in  the  streamwise  spectrum.  This  feature  is  probably  a  result  of  the  low  Reynolds 
numbers  of  the  three  flows.  The  high  wave-niunber  range  of  the  streamwise  spec¬ 
trum  is  associated  with  small  structures  and  has  a  —5  slope  in  the  planar  case.  In 
the  transversely  curved  flows  the  axial  spectra  exhibit  steeper  slopes  in  the  high 
wave-number  range,  indicating  a  weakening  of  structures  in  the  buffer  layer  (Blake 
[1986])  that  contribute  to  the  wall  pressure  fluctuations. 
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In  a  transversely  curved  boundary  layer,  the  volume  of  turbulent  flow  that  must 
be  supported  per  unit  wall  surface  area  is  larger  than  in  the  plane  chzuinel  case  by 
a  factor  of  1  +  S/[2a).  This  is  a  natural  curvature  parameter.  It  was  found  that 
when  scaled  with  the  curvature  dependent  length  scale  given  by 


the  axial  (streamwise)  spectrum  of  the  wall  pressure  fluctuations  of  the  plane  chan¬ 
nel  flow  and  of  the  transversely  curved  flows  collapse  in  the  high  wave  number 
range,  as  shown  in  Figure  4.6.  This  length  scale  is  related  to  the  velocity  scale  de¬ 
fined  in  Section  3.2.3  (Equation  3.2.5,  see  Section  4.6).  Spectra  normalized  with 
£ire  shown  in  Figure  4.7. 

The  two-dimensional  spatial  wall  pressure  spectra  are  shown  in  Figure  4.8.  As  in 
the  planar  case  (see  Fig.  12  of  Choi  &  Moin  [1990])  the  transversely  curved  spectra 
2ire  elongated  in  the  the  sp2inwise  direction.  However,  as  the  curvature  increases  the 
azimuthal  (spanwise)  elongation  of  the  spectrum  increases.  In  addition  the  energy 
of  the  large  structures  (small  Skz  or  6k0)  decreases  with  increasing  curvature. 

4.4  Two-point  correlations 

In  the  two  previous  studies  of  the  weill  pressure  fluctuations  in  boundary  layers 
with  transverse  curvature  (Willmarth  &  Yang  [1970]  and  Willmarth  et  al.  [1976]) 
it  was  shown  that,  as  the  curvature  increases,  the  wall  pressure  becomes  better 
correlated  around  the  cylinder.  This  trend  is  also  evident  in  the  present  flows  as 
shown  in  Figure  4.9  o.  Note,  however,  that  in  viscous  imits  the  azimuthal  (spanwise) 
correlation  length  decresises  with  increasing  curvature  (Figure  4.9  b).  In  the  axial 
two-point  correlations  (Figure  4.10),  the  zero-crossing  point,  Rpp(z/S*)  =  0,  is 
reached  at  z/S*  =  3.9  in  the  6/a  =  5  flow  and  at  z/S*  =  4.6  in  the  S/a  =  11 
flow.  Both  values  are  higher  than  the  value  of  2  reported  by  Willmarth  &  Yang 
(1970)  and  Wilmarth  ei  al.  (1976),  however,  when  compared  with  the  value  of  3.4 
reported  by  Choi  &  Moin  (1990)  for  the  plane  channel,  a  trend  towards  larger  axial 
length  scales  with  increasing  curvature  is  discernible. 
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Figure  4.8  Spectra  of  the  wall  pressure  fluctuations  normalized  by  Ur  and  6 
(8Ak0  =  6 /a,  and  the  scale  of  the  ordinate  is  twice  that  of  the  abscissa);  (a) 
o/a  =  5,  (b)  6/a  =  11.  The  contour  levels  are  logarithmically  distributed  from  10”® 
to  10“^  with  exponent  increments  of  0.3. 
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Figure  4.9  Azimuthal  ^spanwise)  two-point  correlations  of  the  wall  pressure  fluctu¬ 
ations  as  a  function  of  (a)  6  zoid  of  (5)  6a^: - Sfa  =  5  (o'*’  «  43); - 6/a  =  11 

(a+  w  21). 
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Figure  4.10  Axial  (streamwise)  two-point  correlations  of  the  wall  pressure  fluctu¬ 
ations:  - Sfa  =  5  and - Sfa  —  11. 

For  small  separations,  the  wall  pressure  correlation  contours  in  the  planar  case 
(Willmarth  et  al.  [1962],  Choi  &  Moin  [1990])  are  nearly  circulzu:  and,  as  the  sepa¬ 
ration  increases,  the  ratio  of  the  spanwise  to  the  streamwise  length  scales  increases. 
In  contrast  with  the  planar  wall  pressure  correlation  contours,  in  a  transversely 
curved  flow  Willmarth  &  Yang  (1970)  (8/a  =  2  and  a"*"  «  4500)  report  wall  pressure 
correlation  contours  which,  for  large  separations,  are  compressed  in  the  spanwise 
direction.  However,  the  S/a  =  4  measurements  of  Willmarth  et  al.  (1976)  seem  to 
be  inconsistent  with  this  apparent  trend.  Figure  4.11  shows  the  correlation  contours 
of  the  wall  pressure  fluctuations  for  the  plane  channel  flow.  For  large  separation  the 
correlation  contours  are  elongated  in  the  spanwise  direction.  In  contrast,  Figure 
4.12  shows  that,  as  the  curvature  increases,  there  is  a  significant  increase  in  the 
ratio  of  the  axial  to  azimuthal  length  scales  for  all  separations. 
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Figure  4.11  Contour  plot  of  two- point  correlationp  the  wall  pressure  fluctuations 
for  plane  channel  (Choi  &  Moin  [1990])  as  a  function  of  the  streamwise  (rj)  and 
spanwise  (rj)  separations.  The  contour  levels  are  from  0.1  to  0.9  with  increments 
of  0.1. 

4.5  Green’s  Function  Representation 

Further  insight  into  the  spatial  structure  of  the  the  wall  pressure  fluctuations 
czm  be  obtained  by  representing  the  pressure  field  in  terms  of  the  Green’s  function 
of  the  Poisson  operator. 

From  the  analysis  presented  in  Appendix  B,  in  transversely  curved  flows  the 
Fourier  coefficient  of  the  wall  pressure  fluctuations  are  given  by 


^<1+1  , 

p(r)  =  /  G(r,ro)f{ro)rodro+ 
Ja 


ag{r,a)^ 


-  (a  +  l)g{r,a  +  1) 


r=a 


dp 

dr 


r=a-l-l 


(4.5.1) 


where  ro  is  the  radial  position  of  the  source,  /  denotes  the  Fourier  transform  of 
the  source  of  the  pressure  Poisson  equation  (Equation  4.1.1  o)  and  Q  is  the  Green’s 
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Figure  4.12  Contour  plot  of  cwo-point  correlations  of  the  wall  pressure  fluctuations 
for  (a)  Sfa  =  b  amd  for  (b)  6 /a  =  11  as  a  function  of  the  axial  (r+  -  streamwise) 
and  aizimuthal  {a^d  -  spanwise)  sepzirations.  The  contour  levels  axe  from  0.1  to  0.9 
with  increments  of  0.1. 

function  of  the  Laplacian  operator  in  cylindriceil  coordinates  with  derivative  bound¬ 
ary  conditions.  The  global  nature  of  the  pressure  field  is  clear  in  Equation  4.5.1 
due  tc  the  volume  integral  of  G{r,ro)f{ro)  over  the  whole  flowfield.  From  Equation 
4.5.1  the  contributions  to  the  wall  pressure  fluctuations  (p^-)  from  sources  located 
in  the  volume  close  to  the  cylinder  surface  at  a  <  Tq  <  are  given  by 
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Figure  4.13  Fractional  contributions  to  the  root-mean-square  wall  pressure  fluc¬ 
tuations  from  somces  located  at  yo  6  [0,ys]: - plane  channel;  cylinders  with 

- 8 /a  =  5  and - 8 /a  =  11  (j/s  =  —  a).  The  result  is  normalized  by  the 

total  weill  pressure  intensity. 


Pr-  =  /  Q{<i,ro)f{ro)rodro+ 
*  Ja 


aG{a,a)^ 


r=a 


r,Q(o,r,)^ 


(4.5.2) 


r=r. 


Likewise,  the  contributions  of  the  outer  part  of  the  volume  of  the  flow  (r^  <  r©  < 
a  -f- 1)  to  the  wall  pressure  fluctuations  (p_+)  are  given  by 


ra+l  ^ 

Pr+  =  /  0{a,  ro)f{ro)rodro+ 

"  Jrs 


A/  \dp 
rsG{a,rs)-~ 


r=rs 


~(a+  l)G(a,a-h  1)^ 
dr 


(4.5.3) 


r=a+l 
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Figure  4.14  Contour  plot  of  two-point  correlations  of  the  contributions  to  the  wall 
pressure  fluctuations  in  the  plane  channel  from  flow  perturbations  in  the  interval 
(a)  [0,  j/s]  and  (6)  [j/s,!]  (1  denotes  the  chemnel  centerline)  as  a  function  of  the 
streamwise  (rj)  and  spanwise  (rj)  separations.  The  contour  levels  are  from  0.1  to 
0.9  with  increments  of  0.1  and  y*  «  0.2  (yj  w  36). 
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Figure  4.15  Contour  plot  of  two-point  correlations  of  the  contributions  to  the  wal 
pressure  fluctuations  for  ^/a  =  5  from  flow  perturbations  in  the  interval  (a)  [a, 
and  (6)  [rs,a  -f  1]  as  a  function  of  the  eixial  -  streamwise)  and  azimuthal  (a‘*’6 
-  spanwise)  separations.  The  contour  levels  are  from  0.1  to  0.9  with  increments  of 
0.1  and  ys  w  0.2  (yj  =  rj  —  o'*"  w  46). 
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Figure  4.16  Contour  plot  of  two-point  correlations  of  the  contributions  to  the  wall 
pressure  fluctuations  for  ^/a  =  11  from  flow  perturbations  in  the  interval  (a)  [a,  rgl 
smd  (6)  [rs,a  -f- 1]  as  a  function  of  the  axial  (rj  -  streamwise)  and  azimuthal 
-  spanwise)  separations.  The  contour  levels  are  from  0.1  to  0.9  with  increments  of 
0.1  amd  j/s  «  0.2  (yj  =  rj  —  a"*"  a  49). 

Similar  expressions  can  be  obtained  for  the  plane  channel.  The  root-mean-square 
of  the  contributions  to  the  wall  pressure  fluctuations  normalized  by  the  total  wall 
r.m.s.  pressure  is  shown  in  Figure  4.13  as  a  function  of  the  cutoff  rg.  For  the 
two  transversely  curved  flows  as  well  as  the  plane  channel  more  than  80%  of  the 
wall  r.m.s.  pressure  fluctuations  are  produced  by  the  inner  part  of  the  domain 
(0  <  To  <  y*  =  Ts  —  a  a  0.2).  Note  that,  as  the  curvature  increases,  there  is  a  slight 
increase  in  the  wall  pressure  fractional  contributions  from  the  near-wall  (yg  <  0.2) 
part  of  the  flow. 
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Small  separation  correlation  contours  result  from  small  scale  pressure  fluctua¬ 
tions  typicaJly  associated  with  the  inner  part  of  the  boundary  layer,  while  the  large 
separation  correlation  contours  have  significant  contributions  from  large  scale  mo¬ 
tions  in  the  outer  part  of  the  boundary  layer.  This  is  emphasized  in  Figure  4.14 
which  shows  the  plane  channel  two-point  correlations  of  p  -  and  p  +  for  ys  ~  0.2. 
It  is  clear  that  the  spanwise  elongation  of  the  planar  wall  pressure  correlation  con¬ 
tours  is  due  to  disturbances  in  the  outer  part  of  the  layer  (Figure  4.14  h).  Note 
also  that  for  this  cutoff  (j/g  w  0.2)  the  shape  of  the  planar  wall  pressure  correlation 
contours  due  to  the  inner  part  of  the  flow  are  slightly  elongated  in  the  streamwise 
direction  (Figure  4.14  a). 

The  two-point  correlations  of  p^-  and  Pj.+  for  ys  =  Vs  —  a  ^  0.2  are  shown  in 
Figures  4.15  and  4.16  for  8/a  =  5  and  for  8/a  =  11,  respectively.  In  both  cases 
the  two  point  correlations  of  the  wall  pressure  fluctuations  due  to  contributions 
from  the  volume  close  to  the  cylinder  have  a  streamwise  aspect  ratio  similar  to 
that  shown  in  Figure  4.12  for  the  contributions  of  the  whole  flow.  On  the  other 
hand,  in  the  8/a  =  b  flow,  the  shape  of  the  correlation  contours  of  the  wall  pressure 
fluctuations  due  to  the  outer  part  of  the  flow  (t/o  6  [ys  ^  0.2, 1.0],  see  Figure 
4.15  6)  have  nearly  circular  shapes  for  small  separations  and  a  spanwise  orientation 
for  large  separations.  The  same  trend  is  also  evident  in  the  8/a  =  11  flow,  although, 
a  streamwise  orientation  is  discernible  in  the  wedl  pressure  correlation  contours  due 
to  sources  in  the  outer  part  of  the  flow  (Figure  4.16  b).  Wall  pressure  correlation 
contours  due  to  the  outer  part  of  the  flow,  obtained  for  larger  values  of  the  cutoff 
iVs  >  0.2,  not  shown),  become  increasingly  oriented  in  the  speinwise  direction  in 
both  transversely  curved  flows  (8 /a  =  5  and  8/a  =  11),  however,  the  spanwise 
stretching  effect  of  the  outer  part  of  the  flow  on  the  wall  pressure  correlations 
decreases  as  the  curvature  increases. 

The  fact  that,  in  the  three  flows,  the  wall  pressure  correlations  due  to  the  near¬ 
wall  part  of  the  flow  are  elongated  in  the  streamwise  direction  implies  that  the 
near-wall  structures  axe  responsible  for  this  shape.  The  curvatmre  trend  shown  in 
Figures  4.14  through  4.16  is  consistent  with  the  length  scale  parameter  L*  based  on 
Vz  (see  Section  3.5,  Figure  3.25)  which  suggests  that  relative  to  the  plane  channel, 
as  the  curvature  increases,  the  ratio  of  the  streamwise  to  spanwise  length  scales 
increases  near  the  wall  for  both  flows  and  away  from  the  wsdl  in  only  the  8/a  =  11 
case.  Thus,  it  is  conjectured  that,  as  the  curvature  increases,  the  inner  part  of  the 
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flow  plays  an  increasingly  importsint  role  in  determining  the  length  scales  of  the 
wall  pressure  fluctuations. 

4.6  Temporal  Spectra 

The  temporal  spectra  of  the  wedl  pressure  fluctuations  of  the  two  transversely 
curved  turbulent  flows  is  compared  to  their  planau-  counterpart  in  Figure  4.17  a.  As 
expected  from  the  streamwise  one-dimensional  spectrum,  the  temporal  spectrum 
decreases  at  all  frequencies  with  increasing  curvature. 

In  the  wall  pressure  temporal  spectra  measurements  of  Willmarth  &  Yang  (1970) 
and  Willmarth  et  al.  (1976),  a  reduction  in  the  intensity  of  the  high  frequency 
rEinge  was  also  observed.  However,  the  effect  of  curvature  on  the  inner  part  of  the 
boundary  layer  in  their  flow  regime  (high  a'*'),  and  therefore  on  the  high  frequency 
range  of  their  wall  pressure  temporal  spectra  should  have  been  negligible.  This 
apparent  inconsistency  in  the  experimental  data  may  be  related  to  probe  spatial 
resolution  problems,  which  were  not  corrected  for. 

Using  the  curvature  dependent  length  scale  defined  in  Section  4.3  (Equation 
4.3.1),  a  curvature  dependent  time  scale  (S^/ur)  can  be  defined: 


(4.6.1) 


When  scaled  with  this  curvature  dependent  time  scale  the  temporal  spectrum  of  the 
wall  pressure  fluctuations  of  the  two  transversely  curved  flows  collapse  with  their 
pleinar  counterpart  in  the  high  frequency  range  (Figure  4.17  6).  Also  note  that  this 
time  scale  is  related  to  the  velocity  scale  u  defined  in  Section  3.2.3  (Equation  3.2.5) 
by 


6 


(4.6.2) 
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4.7  2-D  Spectra  and  Correlations 

The  frequency /axial  (streamwise)  wave-number  power  spectra,  ^pp(kz.,uj),  for 
the  two  transversely  curved  flows  are  shown  in  Figure  4.18.  As  has  been  the  case 
throughout  this  study  the  intensity  decreases  as  the  curvature  increases  for  all 
frequencies  and  axial  length  sceiles.  In  both  cases  the  narrow  aspect  ratios  of  the 
iso-contours  denote  well  defined  convection  velocities. 

The  space-time  correlation  of  the  wall  pressure  fluctuations  is  given  by  the 
Fourier  transform  of  the  frequency  /  streamwise  wave-number  spectrum, 

^  (4.7.1) 

jfcz  W 

and  is  shown  in  Figure  4.19  for  the  two  transversely  curved  flows.  When  measured  in 
viscous  units,  the  axial  length  scales  of  the  auto-correlation  contours  are  virtually 
unaffected  by  curvature,  however,  the  temporal  extent  increases  with  curvature. 
Again,  clearly  defined  convection  velocities  axe  evident  from  the  oblong  shape  of 
the  contours. 

4.8  Convection  Velocity  and  Taylor’s  Hypothesis 

The  idea  that  structures  in  the  flow  are  convected  at  a  velocity  close  to  the  free 
streeim  velocity  is  an  attractive  one  which  has  received  substantial  experimental 
support  in  planar  boundary  layer  flows.  Experimenteilly,  it  is  typically  easier  to 
measure  temporid  spectra,  from  which  streamwise  spectra  are  obtained  by  invoking 
the  concept  of  a  convection  velocity  and  Taylor’s  hypothesis. 

Several  definitions  for  the  convection  velocity  have  been  proposed  in  the  litera¬ 
ture.  The  most  common  ^l^e  obtained  from  the  frequency  streamwise  wave-number 
spectra,  $pp(fcz,u)),  or  from  its  Fourier  transform,  Rpp{rz,rt). 

For  example,  the  convection  velocity  as  a  function  of  the  axial  (streamwise) 
separation  (r^)  for  a  given  time  delay  is  defined  as  the  ratio  rz/rt^  for  which 
Rppi^zi  is  a  meiximum, 


Vc(r,)  = 


dRpp 

dvz 


n-rtc 


=  0. 


(4.8.1  a) 
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to  0.9  with  increments 


Likewise,  the  convection  velocity  can  be  defined  as  a  function  of  the  temporal 
separation  {rt )  according  to 


Ucirt)  = 


!is 

rt 


dR 


pp 


drt 


{rz,rt) 


=  0. 


(4.8.1  h) 


The  convection  velocities  of  the  wedl  pressure  fluctuations  computed  from  Equa¬ 
tions  4.8.1  are  shown  in  Figure  4.20.  The  convection  velocities  are  lower  than  in  the 
plane  channel.  As  the  axial  separation  increases,  the  convection  velocity  increases 
to  about  O.lUoo  for  6/a  =  5  and  to  about  0.65Z7oo  for  S/a  =  11.  Note,  however, 
that  as  a  function  of  the  temporal  separation  the  convection  velocity  is  practically 
constant  in  both  flows  (Uc  —  0.6£/oo)- 

The  convection  velocity  can  also  be  defined  as  a  function  of  the  streamwise 
wave-number  (Wills  ([1970]), 


Uc(k,}  = 


pp 


du} 


=  0, 


(4.8.2  o) 


U>=Uc 


and  as  a  function  of  the  frequency. 


UJ 


UcH  =  -J-, 


pp 


dk. 


{kz,(^) 


=  0. 


ib  -k 

Kz  — 


(4.8.2  h) 


When  expressed  as  a  function  of  the  axial  wave-number  (Figure  4.21)  the  con¬ 
vection  velocity  is  about  0.617oo  for  large  wave-numbers  in  zdl  three  flows.  Likewise, 
as  a  function  of  frequency  the  convection  velocity  (Figure  4.22)  is  about  0.6517oo 
in  the  curved  flows,  which  is  lower  than  the  value  of  about  O.SUoo  for  the  plane 
channel  (Choi  &  Moin  [1990]). 

In  the  cylinder  flows,  the  various  convection  velocities  show  little  variation  among 
themselves.  In  the  following  a  consteint  value  of  the  convection  velocity,  0.65I7oo, 
is  used  to  sczde  the  temporal  spectrum  (^pp(u>))  into  the  axial  spectrum  {Epp{kz)) 
according  to  Taylor’s  hypothesis: 


^PPci^zc) 


<f>ppi^) 

Uc 


(4.8.3) 
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Figure  4.20  Convection  velocity  normalized  by  Uoo  as  a  function  of  the  (a)  axial 
(streamwise)  separation  {r,/6’)  and  of  the  (6)  temporal  separation  (r/");  cylinders: 
0,  for  6/a  =  5  and  a  ,  for  8/ a  =  11. 
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Figure  4.22  Convection  velocity  normalized  by  Uoo  as  a  function  of  frequency: 
cylinders;  o,  for  6/a  =  5  and  £»  ,  for  6/a  =  11. 


no 


hypothesis  {Uc  —  0.65C/oo) 


Comparison  of  (f)pp{uj)  with  EppJJzzc)  (Figure  4.23)  shows  that  Taylor’s  hypothesis 
is  a  better  assumption  for  low  frequencies,  as  was  also  observed  in  the  plaiiar  case 
(Choi  &  Moin  [1990]). 

The  decomposition  of  the  two-dimensional  spectra  ^pp(kz,u)  and  $pp(fcg,u;)  into 
their  streamwise  and  spanwise  similarity  functions  Fz{kzUc/u})  and  Fe(keUc/<jij),  are 
given  by 


^pp(kz,iu)  =  <l>pp(u:)^^Fz{kzUcH, 

(jO 

^pp(J^9i^)  —  ^pp{^)  ^ — ^~Fg{keUc/^), 


(4.8.4) 


respectively  (Corcos  [1964]).  Unlike  the  plane  channel  where  a  self  similar  behavior 
in  the  spanwise  direction  was  observed,  no  self  similar  behavior  is  apparent  in  either 
Fz(kzUc/i^)  or  FglksUclijj),  as  shown  in  Figures  4.24  and  4.25. 
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CHAPTER  5 


Summary  and  Conclusions 


The  main  objective  of  this  study  was  to  investigate  the  space-time  characteristics 
of  turbulent  Ifows  over  long  towed  cylinders.  The  main  conclusions  are  recapitulated 
in  this  chapter. 

For  a  suflSciently  large  ratio  of  the  cylinder  radius  to  the  boimdary  layer  thickness 
{JSja)  curvature  affects  the  outer  part  of  the  layer.  If  the  iimer  flow  is  to  be  influ¬ 
enced  by  the  cmwatme  then  the  ratio  of  the  cylinder  radius  to  the  viscous  length 
scale  a/(u/ur)  must  also  be  small.  Two  simulations  of  transversely  curved  turbu¬ 
lent  flows  were  successfully  performed.  The  curvature  parameters  were  S/a  =  5 
(a'*‘  «  43)  and  S/a  =  11  (a"*"  «  21). 

In  agreement  with  experimental  measurements,  as  the  curvature  increases,  the 
skin  friction  coefficient  (Cf)  increases  and  the  slope  of  the  mean  velocity  profile  in 
the  logarithmic  region  decreases.  When  a"*"  is  sufficiently  small,  the  viscous  region 
of  the  velocity  profile  deviates  slightly  from  the  planar  law  of  the  wall. 

Measurements  of  the  turbulence  intensities  and  the  Reynolds  shear  stress  (when 
scaled  with  wsJl  variables)  show  a  decrease  with  increasing  curvature  in  the  outer 
part  of  the  flow.  These  outer  flow  characteristics  are  reproduced  by  the  present 
computations.  In  addition,  and  unlike  their  experimental  counterparts,  close  to  the 
wall,  the  computed  turbulence  intensities  and  the  Reynolds  sheeur  stress  decrease 
with  increasing  curvature.  The  reason  for  the  difference  is  attributed  to  the  large 
values  of  a"*"  in  most  experiments,  where  only  the  outer  layer  is  affected  by  the  cur- 
vatxire.  The  reduction  of  the  turbulence  intensities  with  curvature  is  not  the  same 
among  different  components  and,  as  the  cvurvature  increases,  more  of  the  turbulent 
kinetic  energy  is  in  the  streamwise  velocity  fluctuations.  This  is  a  result  of  the 
lower  intercomponent  energy  transfer  due  to  the  lower  pressure  strain  correlations 
in  the  Reynolds  stress  budget  equations  of  the  transversely  curved  flows. 

A  new  local  velocity  scale,  which  is  a  function  of  both  the  curvature  parameter 
S/a  and  the  distance  to  the  wall  y/S  has  been  obtziined  from  the  mean  streamwise 
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momentum  equation.  When  normadized  with  this  new  velocity  scale  some  turbu¬ 
lence  statistics  of  the  plane  channel  and  the  transversely  curved  flows  collapse  in 
the  outer  part  of  the  flow  (j/'*'  >  30).  Near  the  wall  the  new  velocity  scale  is  unable 
to  account  for  all  the  differences  observed. 

As  the  curvature  increases,  the  correlation  between  the  streamwise  and  wall 
normal  velocity  increases  near  the  wall.  In  the  higher  curvature  case  {8 /a  =  11) 
the  correlation  coefficient  decreases  in  the  outer  layer  indicating  a  tendency  towards 
stabilization.  However,  both  flows  computed  axe  fully  turbulent  and  statistically 
steady. 

All  vorticity  intensities  decrease  with  increasing  curvature  throughout  the  layer, 
with  the  exception  of  the  normal  vorticity  intensity,  which  is  not  affected  near  the 
wall  (!/■*■  <  7).  The  axial  (streamwise)  intensity  suggests  that,  as  the  curvature 
increases,  the  near- wall  streamwise  vortices  are  weaker  but  have  a  radius  similar  to 
that  of  the  plane  channel.  Their  position  away  from  the  wall  is  also  not  affected  by 
the  curvature.  As  the  curvature  increases,  the  vorticity  near  the  wall  tends  to  be 
increasingly  oriented  in  the  direction  normal  to  the  wall. 

Another  effect  of  the  transverse  curvature  is  the  longer  streamwise  length  scale 
of  the  low-speed  streaks.  The  mean  spanwise  spacing  of  the  low-speed  streaks  is 
slightly  less  100  wall  units  in  both  transversely  curved  flows  at  y'*'  «  12.  Thus,  as 
the  circumference  of  the  cylinder  decreases  (as  the  curvature  increases),  there  are 
fewer  low-speed  streaks  aroxmd  the  cylinder. 

Near- wall  internal  she2u:  layers  are  common  in  the  two  transversely  ciurved  flows 
computed.  As  the  shear  layers  lift  off  from  the  surface  of  the  cylinder  they  have 
large  spanwise  length  scsdes  relative  to  the  cylinder  radius. 

The  linear  and  nonlinear  sources  of  pressure  fluctuations  are  strongly  reduced  as 
the  curvature  increstses.  The  sources  aissociated  with  streamwise  vortices  are  still 
the  strongest,  however,  as  the  curvature  increases,  a  new  source  of  pressure  fluctua¬ 
tions  associated  with  the  strong  normal  vorticity  fluctuations  becomes  increasingly 
important. 

As  the  curvature  increeises,  the  axial  (streamwise)  wall  pressure  spectrum  de¬ 
creases  for  all  scales,  whereas  the  azimuthal  (spanwise)  wall  pressure  spectrum 
decreases  only  for  small  wave- numbers  (large  scales).  A  curvature  dependent  outer 
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length  scale,  which  increases  with  curvature,  is  proposed  from  geometrical  argu¬ 
ments.  This  length  scale  collapses  the  streamwise  spectra  of  the  wall  pressure 
fluctuations  of  the  two  transversely  curved  flows  studied  with  that  of  the  plane 
channel  in  the  high  wave  number  range. 

Even  though  the  wall  pressure  fluctuations  become  increasingly  better  correlated 
around  the  cylinder,  the  azimuthal  correlation  length,  when  measured  in  wall  imits, 
decreases  as  the  curvature  increases.  On  the  other  hand,  the  axial  (streamwise)  cor¬ 
relation  length  increases  with  curvature.  Unlike  in  the  planar  case,  the  wall  pressiure 
iso-correlation  contours  are  elongated  in  the  streamwise  direction,  for  both  large 
and  small  separations.  Since  the  streamwise  elongation  of  turbulence  structures 
with  curvature  is  most  pronounced  near  the  wall,  it  is  conjectured  that  the  near 
wall  fluctuations  are  more  important  in  determining  the  length  scales  of  the  wall 
pressure  fluctuations.  Fractional  contributions  of  the  flow  (from  inner  and  outer 
layers)  to  the  wall  pressure  intensity  show  a  minimal  effect  of  curvature  relative  to 
the  plane  channel. 

Like  the  axial  (streamwise)  one-dimensional  spectra,  the  temporal  spectra  of 
the  wall  pressure  fluctuations  of  the  transversely  curved  flows  also  decreases  as 
the  curvature  increases.  A  new  curvature  dependent  time  scale  which  increases 
with  increasing  curvature  was  propsed.  When  scaled  with  the  mean  wall  shear  and 
this  time  scale  the  temporal  spectra  of  the  wall  pressure  fluctuations  of  the  two 
transversely  curved  flows  studied  and  that  of  the  plane  channel  collapse  in  the  high 
frequency  range. 

The  two-dimensional  spectra  and  space  time  correlations  of  the  weill  pressure 
fluctuations  give  a  lower  convection  velocity  {Uc  ~  O.GUoo)  than  in  the  plane  chan¬ 
nel.  With  this  convection  velocity  Taylor’s  hypothesis  holds  in  the  two  transversely 
curved  flows  studied. 
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APPENDIX  A 


Reynolds  Stress  Budgets 


In  this  appendix  the  transport  equations  for  the  Reynolds  stresses  and  turbulent 
kinetic  energy  are  presented.  They  are  derived  from  the  Navier-Stokes  equations 
in  general  tensor  notation  (see  for  example  Moser  &  Moin  [1984]) 


I 

I 


ay* 

dt 


+  v‘v<j,  =  -/‘Pt  + 


.Hi 


i  =  1,2,3, 


(A.1) 


where  the  superscripts  denote  contravariant  tensors  and  the  subscripts  following 
a  comma  denote  covariant  derivatives.  The  contravariant  tensor  is  the  metric 
tensor,  which  in  cylindrical  coordinates  is 

..  /I  0  0\ 

-  0  l/r?  0  ,  (A.2) 

\0  0  ij 


and  the  coordinates  (xi,a;2,a:3)  correspond  to  {r,0,z),  respectively. 
The  evolution  equation  for  the  mean  velocity  is  given  by 


.  T7fcT7.- 

dt 


and  the  velocity  fluctuations  equations. 


are  obtained  by  subtracting  Equation  A. 3  from  Equation  i4.1. 


{A.3) 


^  =  -g\i,  +  -  (vV  -  (A.4) 
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From  Equation  >1.4,  evolution  equations  for  the  Reynolds  stresses  are: 


dv^vJ 


^  +V  (VvJ)ji.= 


total  rate  of  change, 
production  rate. 


—  it  +  g^^v*p  jt  velocity  pressure-gradient  term, 


,k 


turbulent  transport  rate. 


viscous  diffusion  rate. 


viscous  dissipation  rate. 


The  transport  equation  for  the  normal  stress,  vl, 


dvl  -  dp 
— t  =  -2vr— 
dt  '^dr 


Her  r  OT  \  or  /  r* 


2 

Rtr  \dr  J  r^\d9  )  \dz ) 


(A.r>) 


(.4.6) 


is  obtained  from  Equation  A.5  by  setting  i  =  j  =  1,  evaluating  the  contravariant 
derivatives  smd  taking  advantage  of  the  homogeneity  in  the  0  and  z  directions 


The  transport  equation  for  the  aaimuthal  stress, 


dt  r  86 


is  similarly  obtained  from  Equation  -4.5  by  setting  i  =  j  =  2. 
The  transport  equation  for  the  axial  stress,  t;|. 


(A7) 


is  similarly  obtained  from  Equation  4.5  by  setting  i  =  j  =  3.  The  budgets  of 
ul,  and  v\  for  B/a  =  5  and  Sfa  =  11  are  shown  in  Figures  A.l,  A. 2  and  A. 3, 
respectively.  The  various  terms  of  the  budget  equations  are  qualitatively  like  their 
planar  counterparts  (Mansour  tt  al.  [1988j)  but  have  their  magnitude  decreased  as 
the  curvature  increases. 
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The  transport  equation  for  the  Reynolds  shear  stress,  VzVr, 


is  similarly  obtained  from  Equation  A.5  by  setting  i  =  1  and  j  =  3.  The  budget 
of  for  6 /a  =  5  and  6 /a  =  11  is  shown  in  Figure  A.4.  As  in  the  case  of  the 
normal  stresses,  the  budget  of  the  Reynolds  stress  tJJiv  is  also  similar  to  its  planar 
counterpeirt  (Mansour  et  al.  [1988])  but  the  magnitude  of  the  terms  decreases  as 
the  curvature  increases. 
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The  transport  equation  for  the  turbulent  kinetic  energy,  =  {\/2){v%-\-v]  +  v%). 


is  obtained  from  the  trace  Ekjuation  A.5. 

The  budget  of  turbulent  kinetic  energy  for  SJa  =  5  and  8/ a  =  11  is  shown  in 
Figure  A. 4.  The  five  terms  in  the  right  hand  side  of  the  transport  equation  for 
the  turbvilent  kinetic  energy  (Ekjuation  A.  10)  of  the  two  transversely  curved  flows 
are  compared  to  their  planar  counterparts  (Mansour  ei  al.  [1988])  in  Figures  A. 6 
through  A.  10.  The  magnitudes  of  ^dl  terms  decreases  as  the  curvature  increases. 
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Figure  A.l  Transport  equation  balance  of  the  eodal  turbulence  intensity  normalized 

by  u  and  Ur  for  (a)  Sfa  =  b  and  (b)  S/a  =  11  :  - production, - dissipation 

rate, . velocity-pressure  gradient  (pressure  diffusion), - turbulent  tremsport, 

- viscous  diffusion. 
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Figure  A. 2  Transport  equation  balance  of  the  normal  turbulence  intensity  nor¬ 
malized  by  V  and  Uj  for  {a)  <J/a  =  5  and  (6)  Sja  —  11  :  - production, - dis¬ 
sipation  rate, .  velocity-pressure  gradient  (pressure  diffusion), - turbulent 

transport, - viscous  diffusion. 


Figure  A. 3  Transport  equation  balance  of  the  azimuthal  turbulence  intensity  nor¬ 
malized  by  V  and  for  [a)  6/a  =  5  and  (6)  6/a  =  11  :  - production, - dis¬ 
sipation  rate, .  velocity-pressure  gradient  (pressure  diffusion), - turbulent 

transport, - viscous  diffusion. 
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Figure  A. 4  Transport  equation  b2kl8mce  of  the  Reynolds  sheaa-  stress  normzdized 

by  V  and  Ur  for  (a)  S/a  =  5  and  (6)  6/a  =  11  ;  - production, - dissipation 

rate, . velocity-pressure  gradient  (pressure  diffusion), - turbulent  transport, 

- viscous  diffusion. 
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Figure  A. 5  Transport  equation  balance  of  the  turbulent  kinetic  energy  normalized 

by  V  and  Ut  for  (a)  6/a  =  5  and  (6)  6/a  =  11  :  - production, - dissipation 

rate, . velocity-pressure  gradient  (pressure  diffusion), - turbulent  transport, 

- viscous  diffusion. 
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Figure  A.6  Velocity  pressure  gradient  term  of  the  turbulent  kinetic  energy  equation 

normalized  by  u  and  Ut: - plane  channel  (Mansour  ei  al.  [1988]);  cylinders  with 

- 6/a  =  5  and - 6/a  =  11. 


Figure  A. 8  Production  of  turbulent  kinetic  energy  normalized  by  v  and  Ut'- 
- plane  channel  (Mansour  et  al.  [1988]);  cylinders  with  -  8/a  =  5  and 


Figure  A. 9  Turbulent  transport  of  turbulent  kinetic  energy  normalized  by  u  and 
Ut-: - plane  channel  (Mansour  ei  al.  [1988]);  cylinders  with - S/a  =  5  and 


Viscous  diffusion 


Figure  A. 10  Viscous  diffusion  in  turbulent  kinetic  energy  equation  normalized  by 

u  and  Ur'. - plane  channel  (Mansour  et  al.  [1988]);  cylinders  with - 6 fa  =  5 

and - S/a  =  11. 
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APPENDIX  B 


The  Green’s  Function  for  Pressure 


In  cylindrical  coordinates  the  Poisson  equation  for  the  pressure  is  given  by  Equa¬ 
tion  4.1.1.  For  each  wave-number  pair  (m,  kz)  the  Fourier  transform  of  Equation 
4.1.1  is  given  by 


<Pp  1  dp 
dr^  r  dr 


1  <Pvr 

dp 

rt  •  1  u*  Vf 

dr 

Rct  dr^ 

r=a  r=a-|-l 

T>  5  ' 

Rer  r^  r 

(B.\) 


where  /(r)  is  the  Fourier  transform  of  the  pressure  source  terms.  In  the  Strum- 
Liouville  form  the  Equation  B.l  is 

^(p)  =  rf(r),  {B.2) 


where  the  operator  C  is  given  by 


(B.S) 


The  Green’s  function,  Q,  of  the  operator  £  with  homogeneous  derivative  bound¬ 
ary  boundary  conditions  is  the  solution  of 


C{Q)  =  Sir  -  ro), 


dr 


r=a 


(B.4) 


r=a-f-l 


=  0, 
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where  6{r  —  Tq)  is  the  Dirac  delta  ’function’  and  Vq  denotes  the  radial  position  of 
the  source.  At  r  =  r©,  ^  is  continuous 


Gir,ro) 


r=r 


+ 

O 


Qir,  ro) 


(B.5) 


and  satisfies  the  jump  condition 


-^^(r,ro) 

r=rt 


r=ro 


ro' 


(B.G) 


For  this  system,  Green’s  identity  is  given  by 


r  [p£(e)  -  e£(p)]  <tr  = 


(.dg  .dp 
''V-d;-^Tr 


a+1 


(B.7) 


which,  upon  integration,  gives  the  general  form  of  the  Fourier  coefficient  for  the 
pressure  fluctuations, 


p(r)=  G{r,ro)f{ro)rodro+ 

Ja 


ag(r,a)^ 


r=a 


-{a  +  l)Q{r,a  +  l)-£ 


r=o+l 


(B.8) 


The  Fourier  coefficient  of  the  wall  pressure  fluctuations  is  obtained  by  evaluating 
this  expression  at  r  =  a. 


i*a+l 

Pw=  Q{a,ro)firo)rodro+ 

Ja 


-  dp 

ag{a,a)£ 


r=a 


-  {a  +  \)g{a,a  +  1)-^ 


r=a+l 


(B.9) 
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The  contributions  to  the  wall  pressure  from  a  voliune  arotmd  the  cylinder  bounded 
by  a  <  To  <  Ts  are  given  by 


Pr-  =  /  Qia,ro)f(ro)rodro+ 

"  Ja 


aQ(a,a) 


dr 


r=a 


-r,S(a,r.)| 


r=r. 


(B.IO) 


Likewise,  the  contributions  of  the  outer  part  of  the  volume  (r^  <  r©  <  a  +  l)  are 
obtained  by  subtracting  Equation  B.IO  from  Equation  B.9 


Jra+l 

'  0(a,ro)f(ro)rodro+ 
rs 


rsQ(a,  ra) 


dp 

dr 


r=rs 


~  (a  +  +  1)^ 


r=a+l 


(B.ll) 
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APPENDIX  C 


Flow  Databases 


Velocity  and  pressure  fields  for  the  two  transversely  curved  flow  calculations 
{S/a  =  5  and  11)  were  saved  every  200  time  steps.  The  time  step  was  2.52  x  lO”"* 
6/ut  time  units  for  the  ^/a  =  5  case  and  2.06  x  10~'*  S/ur  time  units  for  the  6/a  =  5 
case.  The  wall  pressure  and  wall  shear  stresses  were  saved  every  10  time  steps  in 
both  calculations. 

The  data  was  archived  in  Fovuier  space  in  the  axial  and  azimuthal  directions 
and  in  physical  space  in  the  radial  direction  and  in  time.  The  spatial  character¬ 
istics  of  the  data  are  described  in  Table  2.2,  and  the  temporal  characteristics  are 
siunmarized  in  Table  4.1. 

C.l  Velocity  and  Pressure  Data 

The  first  record  of  each  restart  file  is  a  header  containing:  {NR,  NZ,  NTH, 
MTH,  MR,  REY,  MEANPG,  DT,  LTH,  LZ,  RA,  TIME,  TSTEP),  which 
are  the  number  of  planes  in  the  radial  direction  minus  1  {NR),  the  number  of 
Fourier  modes  in  the  axial  direction  {NZ),  the  number  of  Fourier  modes  in  the 
azimuth2d  direction  {NTH),  the  drawer  size  in  the  azimuthal  direction  {MTH), 
the  drawer  size  in  the  radial  direction  {MR),  the  Reynolds  number  based  on  the 
unity  velocity  in  which  the  data  is  reported  (in  these  units  the  friction  velocity  is 
1.2615  for  the  S/a  =  5  case  and  1.3947  in  the  S/a  =  11  case)  and  S/2  {REY)  (in 
both  cases  this  Reynolds  number  is  85),  the  time  step  in  these  units  {DT),  the 
length  of  the  computationeil  domain  in  the  azimuthal  direction  {LTH),  the  length 
of  the  computational  domain  in  the  axial  direction  {LZ),  both  in  units  of  S/2, 
the  ratio  of  the  outer  to  inner  radius  of  the  computational  domain  {RA),  the  time 
{TOLD)  and  the  time  step  {TSTEP)  (in  these  units  the  time  step  was  4  x  10“'*  for 
6/a  =  5  and  3  X  10“'*  for  S/a  =  11).  Next  there  are  {NR/MR-\- 1)  x  {NTH /MTH) 
V D AT A{RGR,THGR\  records  containing  the  complex  Fourier  coefficients  of  the 


137 


velocity  field  in  the  order  (F,,  Vr,  V^)  followed  by  (NR/MR  +  1)  x  (NTH/MTH) 
PDATA[RGR,THGR]  records  containing  the  complex  Fourier  coefficients  of  the 
pressure  field.  The  indices  RGR  and  THGR  identify  the  drawers  of  the  database. 

Each  V  DAT  A  zmd  PDAT  A  drawer  is  a  complex  array  dimensioned  as  DATAV 
[iVZ/2,  MTif,  3,  Mi?]  and  as  DATAP[NZf2,MTH,MR],  respectively.  The  first 
index  indicates  the  axial  (streamwise)  wave  niimber  Ii  =  1,  ...,NZ/2.  The  second 
index  indicates  the  a.zimuthal  (spanwise)  wave  number  —  1,  ...,MTH,  and  the 
wave  number  associated  with  each  is  given  by 

(THGR  -  1 )  ★  MTH  +  h,  h<  MTHf2, 

-[(THGR  -  1)  *  MTH  +  I2-  MTHI2],  h  >  MTH  12. 

The  third  index  in  DATAV  indicates  the  velocity  component  according  to  I3  = 
Vr,3^  V.). 

The  velocity  field  is  given  in  physical  space  in  the  radial  direction  according 
to  the  radial  coordinate  mapping  described  by  Equation  2.1.6.  The  last  index 
indicates  the  radial  plane  RPLANE  =  /<  +  (RGR  —  1)  x  Mi?,  I4  =  1, ...,  Mi?,  and 
$  =  -  cos(ir(RPLANE  -  1)/NR). 

The  restart  files  described  above  are  in  the  mass  storage  system  Columbia,  arc . 
nasa.gov  at  the  NASA  Ames  Rese^ch  Center  in  the  directory  /csf/rf/rft/neves 
/r\m06/rst  for  the  S/a  =  5  calculation  and  in  the  directory  /csf/rf/rft/neves 
/r\ml2/rst  for  the  S/a  =  11  calculation.  The  file  neunes  are  wire.R.  [edititon 
number] . 

C.2  Wall  Data 

There  are  32  files  with  the  wall  data  for  the  6/a  =  5  calculation  and  26  files  for 
the  S/a  =  11  c2dculation.  Each  file  contmns  160  time  samples  of  the  wall  pressure 
(pu,)  and  wall  shear  stresses  (dVi/dr\y,)  and  (dVz/dr\w)-  The  first  record  of  each 
file  contains  (NTI,  HZ,  NTH,  MTH,  MTI,  REY,  LTH,  LZ,  RA),  which  are 
the  number  of  time  saunples  (NT I),  the  number  of  Fourier  modes  in  the  axiad 
direction  (NZ),  the  ntimber  of  Fourier  modes  in  the  eizimuthal  direction  (NTH), 
the  drawer  size  in  the  azimuthal  direction  (MTH),  the  drawer  size  in  the  temporal 
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direction  {MTI)^  the  Reynolds  number  based  on  the  unity  velocity  in  which  the 
data  is  reported  (in  these  units  the  friction  velocity  is  1.2615  for  the  S/a  =  5  case 
and  1.3947  in  the  6/a  =  11  case),  and  6/2  (REY),  the  length  of  the  computational 
domain  in  the  2uzimuthal  direction  (LTH),  the  length  of  the  computational  domain 
in  the  axial  direction  {LZ),  both  based  on  6/2  and  the  ratio  of  the  outer  to  inner 
radius  of  the  computational  domain  (RA). 

Next  there  are  (NTI/MTI)  x  (NTH/MTH)  PDATA[TIGR,THGR]  records 
containing  the  complex  Fourier  coefficients  of  the  wall  pressure  (pw)^  {NTI/MTI)x 
(NTH/MTH)  WDATA[TIGR.,THGR]  records  containing  the  complex  Fourier 
coefficients  of  the  azimuthal  wall  shear  stress  (dVg / dr\y))  and  (NTI/MTI)  x  (NTH 
/MTH)  UDATA[TIGR,THGR]  records  conteiining  the  complex  Fourier  coeffi¬ 
cients  of  the  axial  wall  shear  stress  (914/5r|«;).  The  indices  TIGR  and  THGR  iden¬ 
tify  the  drawers  of  the  database.  Each  PDATA  drawer  is  a  complex  array  dimen¬ 
sioned  as  DATAP[NZ/2,  MTH,  MTI].  The  first  index  indicates  the  cixial  (stream- 
wise)  wave  number  Ii  =  l,,..,NZ/2.  The  second  index  indicates  the  azimuthal 
(spanwise)  wave  number  =  1,  ...,MTH  as  described  in  the  previous  section.  The 
third  index  indicates  the  time  sample  within  the&leTPLANE  =  Is+TIGRxMTI, 
la  =  1,  ...,MTI.  The  time  sequence  through  the  files  is  identified  through  the  file 
edition  number.  The  structure  of  UDATA  and  WDA.TA  is  the  same  as  that  of 
PDATA. 

The  wall  data  files  described  above  are  in  the  mass  storage  system  Coliimbia. 
arc.nasa.gov  at  the  NASA  Ames  Research  Center  in  the  directory  /csf/rf/rft 
/neves/run06  for  the  6/a  =  5  calculation  and  in  the  directory  /csf /rf/rft/neves 
/nml2  for  the  6/a  =  11  calculation.  The  file  names  are  waldat  .R.  [edititon 
number] . 
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